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ABSTRACT 

In  microwave  switching  circuit  design,  a  switching-semiconductor  figure  of  merit 
can  be  defined  and  used  as  the  measure  of  device  switching  effectiveness.  In  this  thesis 
research  switching-semiconductor  figures  of  merit  were  applied  to  general  switching 
circuits.  Computer-aided  circuit  analysis  program  (Touchstone),  and  a  Fortran 
program  were  used  for  the  computation  of  switching  circuit  performance. 
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I.  INTRODUCTION 

A.       BACKGROUND 

The  switching  function  as  performed  in  microwave  circuitry  utilized  in  radar 
applications  is  essential  for  the  performance  of  the  circuit  objective.  In  the  majority  of 
practical  cases,  switching  must  be  done  by  semiconductor  devices,  since  mechanical 
switches  are  ruled  out  by  requirements  of  speed  and  allowable  system  weight  and  size. 
Therefore  the  semiconductor  switch  is  a  vital  component  in  most  microwave  circuit 
applications. 

Despite  its  practical  importance,  relatively  little  research  attention  has  been  given 
to  fundamental  aspects  of  the  switching  problem,  in  comparison  with,  for  example,  the 
problem  of  microwave  amplifier  or  filter  design.  Therefore  it  is  the  purpose  of  this 
thesis  to  examine  some  of  the  basic  features  of  the  switching  problem.  A  general 
requirement  in  the  design  of  a  switching  circuit  is  the  conversion  of  the  set  of  given 
switch  performance  specifications  into  a  set  of  device  requirements  for  the  necessary 
semiconductors  to  be  used  in  the  switch.  It  is  known  that  a  device  figure  of  merit  can 
be  defined  which  provides  a  unifying  factor  in  this  conversion  process.  Some  progress 
in  the  definition  and  use  of  switch  figures  of  merit  has  been  made  by  earlier  workers, 
but  few  practical  procedures  have  been  developed,  and  further  work  is  needed  in  this 
area.  It  is  the  purpose  of  this  thesis  to  investigate  the  characterization  of  several 
aspects  of  switch  performance  as  a  function  of  device  figure  of  merit,  as  well  as  to 
investigate  the  important  aspect  of  the  selection  and  use  of  devices  so  as  to  maximize 
the  power-handling  capacity  of  the  switch. 

The  fundamental  specifications  governing  microwave  switch  selection  and/or 
design  are: 

1.  Low  insertion  loss  in  on-state 

2.  High  isolation  in  off-state 

3.  Broad  bandwidth 

4.  SulTicient  power  capability 

These  characteristics  are  interactive,  a  change  in  one  affecting  the  others.  The 
microwave  switches  considered  here  are  those  which  employ  semiconductor  devices 
(p-i-n  diodes  and  FET's)  operated  in  two  alternate  bias  states.  The  device  is 
characterized  by  a  pair  of  two-terminal  impedance  values  Zj  and  X^  in  the  respective 


bias  states.  The  parasitic  capacitance  associated  with  semiconductor  switching  devices 
necessitates  some  form  of  impedance  compensation  in  order  to  achieve  acceptable 
insertion  losses  for  the  switch.  Consequently  the  theoretical  analysis  of  switching 
circuitry  has  been  oriented  toward  the  study  of  two-state  impedances  embedded  in 
lossless  impedance  transformers  [Refs.  1.2,3].  These  researches  on  switching  theory 
have  introduced  specific  invariant  parameters  related  to  the  device  impedance  in  its  two 
states,  under  lossless  impedance  transformation. 

Kurokawa  and  Schlosser  [Ref.  1],  showed  that  diodes  incorporated  in  a  one-port 
reflective  switch  can  be  characterized  by  a  quality  factor  parameter,  Q"^  : 

q2      =  r    21  (        11) 

R1R2 

where     Ri  is  the  resistive  component  of  Zp  and 

R2  is  the  resistive  component  of  Z2. 
The  numerical  value  of  Q  is  unchanged  when  the  impedances  are  subjected  to  a  lossless 
impedance  transformation. 

Atwater  and  Sudbury  [Ref  2]  showed  that  this  Q  had  the  character  of  a  figure  of 
merit,  serving  as  an  indicator  of  switching  performance  for  semiconductors  embedded 
in  variety  of  circuit  configurations,  and  not  restricted  to  the  original  reflection-type 
circuit  of  Kurokawa  and  Schlosser.  They  followed  the  research  of  Kurokawa  and 
Schlosser  and  assumed  that  Zi  and  Z2  had  been  transformed  to  real  (resistive)  values 
in  both  states.  Equation  1.1  becomes: 

q2     =  ^    ^'     ^^  (eqn  1.2) 

R1R2 

From  equation  1.2,  we  have  the  relationships: 

Ri  is  a  function  of  Q  and  R-).  and 

R->  is  a  function  ot  Q  and  R, . 
Knowing  the  value  of  Q  and  given  that  Ri,  R2  are  kno'wn,  the  switching  performance 
is  fully  specified. 


Kawakami  [Ref.  3]  found  a  different  form  of  switching  invariant  parameter.  He 
defined  this  parameter  as  a  semiconductor  figure  of  merit,  M: 

IZ1-Z9I 
M      =  ^     -V.  (eqn  1.3) 

IZ^-Z,  1 

The. parameter  M  assumes  values  between  0  and  1,  and  is  also  unaltered  by  lossless 
impedance  transformation.  Kawakami  also  showed  that  when  a  device  was 
incorporated  in  a  symmetrical  two-port  junction  as  a  switch,  its  transmission 
coefficients  in  the  on  and  off  states,  S2][(l)  and  S2i(2),  respectively,  are  related  by 
equation: 

|S-,,(1)1(1-M) 

|Soi{2)l     =  ^^        ■ (eqn  1.4) 

21^  l-M[|l-2S2i(l)|] 

Equation  1.4  fully  characterizes  the  switching  performance  of  the  device,  uith 
parameter  M  as  a  figure  of  merit.  It  may  inverted  to  yield  the  form: 

|S2i(2)l-|S2i(l)l  ,        ,  ,. 

M        =  — — (eqn  1.5) 

|S2i(2)|  +  |S2i(l)|[l-21S2i(2)|] 

This  expression  is  based  on  the  assumption  that  the  switching  device  in  its  two  states 
has  the  properties: 

Fj     =     0  (matched;  Z  =  Zq  ) 

^2     =     ^max     =     ^^ 

Equation  1.5  presents  a  useful  tool  for  the  selection  of  a  semiconductor  device  to 
satisfy  low  power  switching  applications.^ 

A  fundamental  zheorem  was  denved  by  Hines  [Ref  4],  to  predict  the  power  that 
can  be  switched  by  a  singie-diode  switch  when  it  is  connected  in  a  iumped-eiement 
circuit  such  as  to  resonate  the  parasitic  reactance  of  the  switch  at  a  single  frequency. 
Hines  showed  that  the  power  which  can  be  switched  by  the  diode,  P^.^,  is  given  by: 


^In  the  absorption  switch  assumed  in  equation  1.4,  the  diode  absorbs  the  incident 
energy  in  the  switch  ofT-state. 
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Psw      =       0.5[V^,,I^,J  (eqnl.6) 


sw  "•^'-   max'max 

where  V^^^^  is  a  diode  breakdown  voltage,  and 

^max  ^^  ^  maximum  forward  current  rating. 
Hines  conjectured  that  P^^^  represents  a  theoretical  limit  for  power  switching  by  a 
given  device  in  an  impedance-matched  circuit. 

In  many  applications,  microwave  switching  circuits  typically  require  multiple 
diodes  in  cascade  connection  to  achieve  adequate  levels  of  off-state  isolation.  For 
example,  radar  and  phased-array  modules  need  a  double-throw  switch  of  high  isolation 
to  protect  the  receiver  front  end. 

B.       SCOPE  OF  THE  THESIS  RESEARCH 

For  the  application  to  microwave  switching,  the  proposed  semiconductor  figures 
of  merit  are  tested  for  different  types  of  switching  circuitry,  to  find  the  best  parameter 
to  use  in  a  switch  selection  and/or  design  problem  in  order  to  minimize  insertion  loss  in 
the  on-state,  and  maximize  isolation  in  the  off-state.  Additional  objectives  are  to 
determine,  for  the  designed  switch,  what  are  the  power  dissipation  and  heat  removal 
problem  in  the  individual  semiconductor  components. 
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II.  CIRCUIT  ANALYSIS 

In  the  present  chapter  we  generalize  the  circuit  representation  of  microwave 
switches  of  our  interest.  Following  this,  a  suitable  circuit  model  may  be  selected  for 
the  particular  type  of  switch  required.  The  simplest  representation  is  that  of  an  N- 
port,  with  scattering  matrix  description.  The  number  of  ports  is  determined  by  the 
number  of  semiconductor  devices  and  the  type  of  switch  (single-throw  or  multiple- 
throw).  For  the  semiconductor  devices,  there  are  three  possible  ways  of  connection  to 
our  circuit  model,  in  shunt  or  in  series  or  a  combination  of  both.  Because  of  heat- 
sinking  considerations  we  select  the  shunt  connection  and  use  a  one-port  network  to 
represent  the  two-terminal  semiconductor  device. 

A.       CLASSIFICATION  OF  MICROWAVE  SWITCHES 

1.  Single  Pole  Single  Throw  Switches  (SPST) 

The  single  pole*  single  throw  switch  is  a  switch  which  has  one  input  port  and 
one  output  port.  In  this  topology  we  consider  our  switch  as  a  'box'  and  it  is  a  two-port 
box  network  shown  in  Figure  2.1(a).  For  an  ideal  switch,  all  the  power  is  transferred  to 
the  load  at  the  output  port  in  its  on-state,  and  no  power  transferred  to  the  load  in  its 
off-state. 

2.  Single  Pole  N  Throw  Switches  (SPNT) 

For  these  types  of  switches,  there  is  an  input  port  and  N  output  ports.  Their 
switching  circuit  models  are  a  three-port  box  network  for  single  pole  double  throw 
switches,  shown  in  Figure  2.1(b),  and  a  (N+l)-port  box  network  for  single  pole  N- 
throw  switches,  shown  in  Figure  2.1(c).  In  microwave  switch  apphcation,  we  are 
interested  in  single  pole  double  throw  switches  as  mentioned  in  the  introductoiy 
chapter.  For  SPDT  switches,  when  one  arm  of  the  switch  is  in  the  on-state,  the  other 
arm  is  in  the  off-state,  and  vice  versa.  In  general,  we  may  consider  a  SPDT  switch  as 
two  SPST  switches  with  dieir  Input  arms  tied  together  and  apply  two  aiilerent  bias 
states  to  the  semiconductor  devices  m  the  diiTerent  arms  of  the  switch. 


^Single  pole  is  a  line-over-ground  circuit. 
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B.       SCATTERING 
SWITCHES 


MATRIX        REPRESENTATION        OF       MICROWAVE 


We  can  use  the  properties  of  the  S-matrix  [Refs.  5,6:  pp.   112,  50],  defined  in 
Figure  2.2,  and  assume  that  our  switches  are  ideal  switches  mth  matched  load  at  both 
input  and  output  ports.  We  can  represent  the  scattering  matrices  for  both  on  and  off 
states  of  our  switches  as  follows. 
I.  Single  Pole  Single  Throw  Switch 

ON  STATE:  [S]      = 


OFF  STATE:   [S]     = 


0 

1 

1 

0 

1 

0 

0 

1 

2.  Single  Pole  Double  Throw  Switch 


ARM  (2)  ON:     [S]    = 


0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

ARM  (3)  ON:    [S]    = 


C.       SWITCHING  CIRCUIT  MORPHOLOGIES 

1.  Number  of  Semiconductor  Devices 

In  an  SPST  switch  or  one  arm  of  SPDT  s^v^itch  there  is  at  least  one 
semiconductor  device.  For  improvement  of  load  isolation,  the  switch  may  require 
additional  semiconductor  devices.  In  general,  a  SPNT  switch  may  use  M 
semiconductor  devices. 


M 


=     pN 


:eqn2.I) 


where    M  is  the  number  of  semiconductor  devices, 
N  is  the  number  of  throws,  and 
p  is  a  positive  integer. 
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2  PORT 
EOX 


a.   SPST 


b.   SPDT 


N+1   POR:! 
BOX 


4-4  4-'^n■l 


■*2 


■*z 


C.   SPNT 
Figure  2.1     Single  Pole  N  Throw  Switches. 


^    ^IL 


portl 


That  is 
b 


N  PORT 
NETWORK 


b2  .^      ^a2 


■? 


V2      port2 


port  n 

Cb]        =        [s]    [a] 

=     Sii^i     +     512^2  +  - 


~  + 


2ln% 


n 


=     S^^a.^     +     S^2^2  + 
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Fisure  2.2    Scattering  Parameters. 
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2.  Switch  ''BOX"  Topologies 

We  can  define  the  switching  circuit  as  a  lossless  embedding  network  with  one 
input  port,  N  output  ports,  and  M  ports  to  which  the  semiconductor  devices  are 
coupled.   The  total  number  of  ports  of  the  embedding  network  is  given  by: 

Number  of  ports      =       (N+l  +  M)  •  (eqn  2.2) 

Example  of  switching  circuits. 

1.  SPST  switch  with  1  semiconductor  device.  Figure  2.3(a) 

Number  of  ports     =     (1  +  1  +  1) 
=      3     ports 

2.  SPST  switch  with  2  semiconductor  devices.  Figure  2.3(b) 

Number  of  ports     =     (1  +  1  +  2) 
=      4     ports 

3.  SPDT  switch  with  2  semiconductor  devices.  Figure  2.3(c) 

Number  of  ports     =     (2+1  +  2) 
=      5     ports 

D.      SWITCHING  SEMICONDUCTOR  DEVICE  CHARACTERIZATION 

The  semiconductor  device  in  shunt  connection  can  be  represented  by  a  one-port 
network,  as  stated  above.'^  The  device  itself  in  the  switching  circuit  may  be 
characterized  by: 

1.  Two-states  of  impedance  [Ref  2].  Z^  (  i  =  1,2  )  or 

2.  Two-states  of  reflection  coefficient  [Ref  1].   F-   (  i  =  1,2  ). 
The  relations  between  impedance  and  reflection  coefficient  are: 

1  +  T; 

Zq  (eqn  2.3) 


(eqn  2.4) 


^i 

l-Ti 

ri 

= 

Z.;-  Z,1 

Zi-Zo 

3 


The  device  is  in  shunt  to  ground  in  the  single  pole  circuit. 
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a.  SPST  1  DEVICE 


•  2 


♦  2 


b.  SPST  2  DEVICES 


c.    SPOT  2  DEVICES 
Figure  2.3     ExamDie  of  Switching  Circuits. 
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Figure  of  merit  of  the  semiconductor  device  is  defined  based  on  the  device 
characterization  chosen. 

1.  Figure  of  Merit,  M 
From  equation  1.3 


M*      = 


IZ1-Z2I 
|Zi-Z2"1 


Kawakami  postulated  that  the  fundamental  characterization  was  in  terms  of 
two  specific  states  of  F-  .   He  assumed  the  following  conditions: 

r^    =    0       (  Z  is  matched;  Z  =  Zq  ) 

A  lossless  matching  network  is  assumed  to  be  used  to  transform  the  device 
impedance  to  the  matched  condition  in  this  state.  In  the  second  bias  state  of  the 
device,  the  reflection  coefficient  presented  by  the  transformed  device  impedance  is  then: 

''"2'    ~    l^max'  *>  ^^^^  ^^  equal  to  the  figure  of  merit,  M,  of  the  devices  ) 

However,  this  type  of  switch  is  a  reflection  switch.'^ 
2.  Quality  Factor,  Q 

Kurokawa  and  Schlosser  postulated  that  the  fundamental  characterization  was 
in  terms  of  two  states  of  impedance  and  their  real  parts  (resistive).  They  defined  a 
quality  factor  Q  ,  shown  in  equation  1.1,  as  the  figure  of  merit  of  the  semiconductor 
devices. 

From  the  parameters  defined  above  we  can  characterize  the  p-i-n  diode  and 
specify  a  diode  quality  factor  or  the  figure  of  merit  for  the  device  used  in  microwave 
switching  circuits.    Diode  quality  factor,  Q,^iode'  ^'^^  ^  typical  switching  diode  is  equal 
to  the  switching  quality  factor,  as  may  be  shown: 
From  equation  1.1 


Q^ 


R^R2 


'^The  F:    =    0   state,  absorbs  energy  and  is  not  suitable  for  high  power  switching 
circuits. 
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For   a   p-i-n   diode,   a   simple   circuit   model  is   assumed  in  which   the   two 
impedance  states  are: 

Z|         =     Rj  forward  bias 

Z2        =     R2  +  iX2  reverse  bias 

If  we  assume  that   R|  =  R2  =  R  (  for  both  bias  states  ) 

R^ 
1 

(  RCco  )^ 
Thus  Kurokawa-Schlosser's  Q  reduces  to  the  conventional  diode  Q: 

Q^      =     Q^diode  (^q^  2.6) 

3.  Comparison  of  Figure  of  Merit 

We  can  find  the  relationships  between  the  figure  of  merit,  M,  and  the  quality 
factor,  Q,  as  follows: 

Kurokawa-Schlosser,  quality  factor,  from  equation  1.1 


Q 


2        ^     IZi-  Z2 
RjR2 


i2 


Kawakami,  figure  of  merit,  from  equation  1.3 

7               iZi-Z^P 
|M|2      =      U 2L^ 

iz^  +  Zt  r 
-)  (Ri  -  R.r  ^  (X,  -x^r 

|Mr      =      -^-^ } 0  (eqn  2.8) 

(Rj  +  R2)^  +  (Xj  -  X2)2 


"MmP  ,      IZj  -  Z2i^ 

IZi  +  Z2  |2 
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= 

IZi  +  z^V 

IZj  4-  Z{\^  •  |Zi  -  Z2P 

= 

|Z|  +  Z2V 
4R^R2 

» 

= 

q2 

4|M|2 

q2 

= 

4^ 

IMp 

= 

Q2 

4  +  q2 

If   q2  >  >  1 

|M| 

=  ■•? 

(eqn  2.9) 


(eqn2.10) 


(eqn  2.11) 


E.       TRANSMISSION   LINE  SECTION  REPRESENTATION   OF  SWITCHING 
CIRCUIT  MODEL 

From  an  N-port  box  representation,  the  semiconductor  devices  are  in  shunt  to 
ground  and  have  two  alternate  impedance  values.  For  practical  microwave  integrated 
circuit  (MIC)  applications,  it  is  convenient  to  use  transmission  line  sections  in  the 
embedding  network  of  the  switches 

Figure  2.4,  shows  the  circuit  model  of  SPST  and  SPDT  switches  with  one 
semiconductor  device  in  each  arm.  To  optimize  the  switch  characteristics  for  low 
insertion  loss  and  high  isolation,  the  impedance  values  of  the  devices  are  kept 
unchanged  at  rhs  particular  frequency  in  each  state.  The  length  and  characteristic 
impedance  of  the  transmission  ime  section  are  vaned  to  achieve  the  optimum  switching 
characteristics  in  a  particular  frequency  interval. 
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0     =  Transmission  Line  Electrical  Lengt±i 
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Fi<'ure  2.4    Transmission  Line  Model  of  Switch  Networks. 
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F.       APPROACHES  TO 
DISSIPATION 


COMPUTATION 


OF        CIRCUIT        POWER- 


Using  the  Touchstone  program,  we  can  design  a  switch  to  have  optimum 
switching  characteristics.  But  the  semiconductor  devices  in  the  optimized  circuit  may 
not  be  able  to  handle  the  amount  of  power  dissipated  due  to  heating  of  the  devices. 
For  effective  switching  design,  it  is  advisable  to  calculate  the  maximum  power 
dissipation  in  each  device  before  construction  of  the  switch. 

1.  ABCD  Matrix  Representation  of  a  Two-port  Network. 

From  Figure  2.5  the  ABCD  matrix  is  defmed  in  equations  2.12  and  2.13 


Vi    =     AV2  +  BI2 

Ij    =     CV2  +  DI2 
Equations  2.12  and  2.13  have  the  matrix  form. 


A      B 
C     D 


^2 


(eqn2.12) 
(eqn2.13) 

(eqn2.14) 


"A 

B" 

-1 

vf 

.C 

D. 

■hi 

'D 

-C" 

■vf 

_-B 

A. 

Uii 

(eqn2.15) 


The  utility  of  the  general  circuit  parameter  matrix  lies  in  the  fact  that  cascaded 
two-ports  can  be  represented  by  the  chained  product  of  element  matrices  [Ref  5:  p 
110]. 


= 

"Ai     b,' 

52      D2_ 



An      B„' 

-In+lJ 

(eqn  2.16) 


Thus  the  general  circuit  parameter  matrix  representation  of  cascaded  two- 
ports  is  the  matrix  product  of  general  circuit  parameter  matnces  of  the  individual  two- 
ports,  as  shown  m  Figure  2.6. 

2.  ABCD  Matrix  Representation  of  the  Switching  Circuit  Model 

In  the  circuit  model  we  have  transmission  line  sections  and  semiconductor 
devices  combined.  Lumped-circuits  can  be  characterized  by  ABCD  two-port  matrices 
[Ref  7:  p.  188],  as  follows: 
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Figure  2.5    ABCD  Matrix  Model. 
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Figure  2.6    N  Two-Port  in  Cascade  with  Equivalent  Matrix. 
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2.1  Series  impedance,  as  shown  in  Figure  2.7  (a) 


A      B 
C      D 


1      Z 
0      1 


(eqn2.17) 


2.2  Shunt  admittance,  as  shown  in  Figure  2.7  (b) 


A      B 
C      D 


1      0 
Y      1 


(eqn2.18) 


2.3  Lossless  transmission  line  section,  as  shown  in  Figure  2.7  (c) 


A      B 
C      D 


CosG  jZQSinG 

jYQSinG      Cose 


(eqn2.19) 


0  is  an  electrical  leneth  of  a  transmission  line  section 


2.4  Lossy  transmission  line  section,  as  shown  in  Figure  2.7  (d) 


A      B 
C     D 


CoshCY  1)         ZoSinh(7  1 ) 
YQSinh(Y  I)     Cosh(Y  I) 


(eqn  2.20) 


1   =    Physical  length  of  a  transmission  line  section. 
Y   =    a  +  jp 

a    =    transmission  line  attenuation  constant. 
P    =    wavenumber. 
3.  Calculation  of  Power  in  a  Cascade  Circuit 
From  Figures  2.5,  and  2.6  [Ref  8:  p.  393]. 


Z: 


AZl  +  B 


m 


CZl  + 

D 

v. 

Zg- 

Zin 

Zin 

-^ 

V 

Z    4-  Z-     g 
^2      ^m 


Ij+1        =     (Ajlj)  -  (CjVp 


(eqn  2.21) 

(eqn  2.22) 

(eqn  2.23) 
(eqn  2.24  ) 


^  a    =  0  for  a  lossless  transmission  line  section. 
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Figure  2.7    ABCD  Matrix  Representation  of  Devices. 
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Vj+1        =    (DjVp  -  (Bjip  (eqn2.25) 

i    =  1,2,3, ,n 

Power  dissipated  in  each  two-port  network. 

?^         =    0.5^e[(Vi^Ii^*)-(V]^+lIk+l*)]  (eqn2.26) 

k   =  1,2,3, ,n 

Power  at  load 

Pl         =    O.S^elY^^yl^^^*]  (eqn2.27) 

=    Pin-f  Pk  ^^^^2.28) 

k=l 
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III.  MICROWAVE  SWITCHING  DIODE 

The  use  of  the  semiconductor  diode  as  a  switching  element  in  microwave  circuits 
is  based  on  the  difference  between  the  diode  impedance  under  reverse  and  forward  bias. 
The  diode  appears  as  a  very  small  impedance  under  forward  bias  and  a  very  large 
impedance  under  reverse  bias.  For  a  shunt-connected  diode  as  a  microwave  switch,  in 
its  on  state  most  of  the  power  is  delivered  to  the  load,  and  the  diode  looks  like  a  large 
impedance  (reverse  bias  state).  On  the  other  hand,  in  its  off  state,  the  power  is 
reflected  back  to  the  source,  or  dissipated  in  the  diodes,  and  the  diodes  look  like  a 
small  impedance  (forward  bias  state)  in  shunt  to  ground.  In  this  chapter  we  will 
introduce  some  characteristics  of  the  p-i-n  diode,  its  behavior  under  forward  and 
reverse  bias  states,  the  equivalent  circuit,  and  power  handling  capability. 

A.       THE  P-I-N  DIODE  STRUCTURE 

In  a  p-i-n  diode,  the  semiconductor  wafer  has  a  heavily  doped  P-region  and  a 
heavily  doped  N-region  separated  by  a  layer  of  high  resistivity  material  that  is  nearly 
intrinsic.  The  thickness  of  the  high  resistivity  layer  lies  in  the  range  between  1.5  to  150 
microns  [Ref  9:  p.  5].  The  electrical  contacts  are  made  to  the  two  heavily  doped 
regions. 

In  practice,  there  are  no  materials  without  any  impurities.  A  practical  p-i-n 
diode  consists  of  extremely  high  resistivity  P  or  N-type  material  between  low  resistivity 
(highly  doped)  P  and  N-regions.  The  lightly  doped,  high  resistivity  layer  is  called  7t-type 
for  P  material,  and  v-type  for  N  material.  The  I-region  of  a  p-i-n  diode  can  consist  of 
either  V  or  TC-type  material.  Figure  3.1(a)  shows  a  P-v-N  diode  structure.  If  the  I- 
region  is  of  sufficiently  high  resistivity,  with  few  impurity  atoms  that  will  be  ionized 
the  depletion  zone  will  extend  throughout  the  I-region  and  include  a  small 
penetration  into  both  the  P  and  N-regions.  Because  of  heavy  doping  in  the  P  and  N- 
regions  [he  depieiion  zone  wiU  not  extend  ven.^  far  :nto  ihem,  ana  che  vvidth  ot  the 
depletion  zone  will  be  essentially  equal  to  the  I-layer  width,  \V.  Figure  3.l{b)  shows 
the  structure  of  P-Ti-N  diode,  the  depletion  zone  width  is  approximately  equal  to  the 
width  of  the  intrinsic  region.  Most  p-i-n  diodes  use  v  material  [Ref  7:  p.  41]  for  the 
I-laver. 
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B.       THE  ELECTRICAL  EQUIVALENT  CIRCUIT  MODEL 

The  p-i-n  diode  used  in  microwave  application  has  a  simplified  equivalent  circuit 
under  forward  and  reverse  bias  states  as  showTi  in  Figure  3.2. 
L  Forward  Bias  Equivalent  Circuit 

Under  the  forward  bias  state,  both  holes  and  electrons  are  injected  into  the  I- 
region,  and  the  mobile  carriers  form  a  conductive  plasma  whose  resistivity  depends  on 
the  carrier  density  obtained  from  the  injection  process.  The  forward  bias  resistance,  Kc 
may  be  expressed  as: 

Rf       =      R^  +  R^^  (eqn3.1) 

where   R-     =    resistance  due  to  the  I  region,  and 

Rj|-  =    resistance  of  P  and  N-regions  and  their  ohmic  contacts. 
For  an  abrupt  junction  diode,  the  resistance  due  to  the   I-region  can  be 
expressed  as  follow: 

R:       =      (eqn  3.2) 

2ifH,,iL 

where    W      =    I-region  thickness, 

i|>     =    DC  forward  bias  current, 

]i^y  =    the  average  carrier  mobility,  and 

Tj^     =    the  carrier  lifetime  in  the  intrinsic  region. 

In  the  shunt  configuration  switching  application  the  forward  bias  resistance, 
R-,  is  the  major  factor  of  power  dissipation  in  the  diode.  Since  R-  is  an  inverse  function 
of  a  DC  bias  current,  i^  ,  we  can  increase  i^  to  minimize  R-  .  For  low  voltage  p-i-n 
diodes,  an  i|-  of  10  mA,  and  for  high  voltage  p-i-n  diodes  150  to  200  mA  is  sufficient 
to  make  R^  equal  to  R^^fRef  9:  p.  8]. 
2.   Reverse  Bias  Equivalent  Circuit 

Under  die  reverse  biased  state,  the  I-region  is  depleted  of  carriers  and  "he  p-i-a 
diode  appears  as  a  constant  capacitance  to  a  microwave  signal.  The  dissipation  losses 
occur  in  the  ohmic  contacts  and  the  resistances  of  the  P  and  N-regions.  R^.  is  used  to 
represent  these  losses  in  the  equivalent  circuit  model.  Due  to  the  capacitive  reactance 
at  microwave  frequency,  the  circuit  impedance  of  the  p-i-n  diode  in  its  reverse  bias 
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Figure  3.1     Profiles  for  the  two  P-I-N  Diodes  Types. 
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Figure  3.2    P-I-N  Diode  Chip  Equivalent  Circuit. 
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state  is  much  larger  than  the  impedance  in  the  forward  bias  state.  The  capacitance,  C: 
of  the  p-i-n  diode  is  determined  by  the  diode  size  and  I-region  thickness  as  shown  in 
equation  3.3 

C:        =      -^^ (eqn3.3) 


where  £q     =    free  space  permittivity, 

Ej.     =    the  relative  dielectric  constant  of  the  material, 
D     =    the  junction  diameter,  and 
W     =    I-  region  thickness. 
This  equation  is  valid  for  all  reverse  bias  voltage  magnitudes  greater  than  the 
minimum  voltage    for   complete    I-layer   depletion.    Figure   3.3    shows   equation   3.3 
graphically  for  typical  available  p-i-n  diode  I-regions. 
3.  P-I-N  Diode  Cutoff  Frequency 

Hines  [Ref  4],  defmed  a  figure  of  merit  for  p-i-n  diodes  to  relate  the  switching 
effectiveness  in  terms  of  switching  cutoff  frequency,  f^^.  The  Hines  figure  of  merit  is 
given  by 


1 

27rC:VRfRj. 


^cs      =     ._^    /^  ^  (eqn3.4) 


The  figure  of  merit,  f^^  is  one  of  the  parameters  used  in  microwave  switching 
design,  and  for  R^-  =  R^.,  f^^g/f  is  equal  to  the  conventional  diode  Q.  Knowing  C:,  Rp 
and  Rj.  we  can  estimate  the  operating  frequency  range  of  the  selected  p-i-n  diode  from 
equation  3.4 

C.       POWER  HANDLING 

In  the  off-state  of  the  microwave  switch,  the  RF  incident  power  is  reflected  back 
to  ihe  scarce  or  dissipated  m  the  p-i-n  diode.  This  dissipation  causes  the  temperature 
of  the  junction  lo  rise.  For  each  diode  there  is  a  maximum  temperature  that  it  can 
withstand.  The  other  factor  which  limits  the  amount  of  power  handling  is  the  peak 
voltage  that  appears  across  the  junction,  which  can  cause  damage  to  the  junction.  In 
the  p-i-n  diode  specification  sheets  these  important  parameters  are  available  for  the 
designer  to  estimate  the  power  handling  capability  of  the  diode  in  the  designed 
switches. 
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1.  Maximum  Voltage  across  the  P-l-N  Diode 

The  instantaneous  voltage  across  a  p-i-n  diode  junction  is  the  algebraic  sum  of 
the  RF  voltage  and  DC  bias  voltage.  If  the  junction  voltage  exceeds  the  device's 
breakdown  voltage,  the  p-i-n  diode  is  in  the  avalanche  region,  and  the  p-i-n  diode  acts 
as  a  low  resistance  causing  high  negative  current  flow  and  rapid  temperamre  rise. 
Figure  3.4  shows  the  condition  for  a  high  RF  voltage  drive  with  reverse  bias.  The 
maximum  voltage  across  the  junction  occurs  when  the  diode  is  in  the  reverse  bias  state. 
The  maximum  instantaneous  voltage  appUed  to  the  junction  is  defined  as: 

Vmax     =      ^RFp,,i,  +  ^bias  (^^^  3.5) 

where  Vj^^^^g  is  the  DC  reverse  bias  impressed  on  the  diode. 

White  [Ref  7:  p.  155],  suggested  the  conservative  limit  of  the  safe  maximum 
voltage  across  the  junction  as: 

Vax     =      0-5  Vfe  (eqn3.6) 

where  V^^  is  the  p-i-n  diode  breakdown  voltage.  A  circuit  designer  may  choose  a 
different  safety  factor  within  an  acceptable  margin  of  safety  for  the  designed  switch. 

2.  Maximum  Junction  Temperature 

The  junction  temperature  of  a  p-i-n  diode  is  determined  by  the  ambient 
temperature  of  the  circuit  and  the  power  dissipated  in  the  diode.  The  dissipated  power 
may  be  pulsed  or  CW.  In  either  case,  junction  temperature  must  be  kept  below  the 
maximum  allowable  value.  This  maximum  temperature  is  dependent  on  the  diode 
material  (semiconductor),  and  the  materials  used  in  mounting  and  connecting  to  the 
semiconductor.  The  recommended  maximum  junction  temperature  for  safe  operation  is 
200 °C  [Ref  10:  p.290]. 

a.    CW  Power  Dissipation 

In  the  case  of  a  continuous  or  a  long-duration  pulse,  ^he  temperature  !imit 

is  determined  by  the  initial  temperature  of  the  diode  matenals,  the  temperature  of  the 

heat  sink,  and  thermal  resistance  between  the  wafer  and  the  heat  sink  [Ref  10:  p. 290]. 

6:„  is  the  total  thermal  resistance  unit  in  °C  /watt.  The  junction  temperature,  T:,  for  a 
J*-  J 

given  power  dissipated,  P^  and  the  ambient  temperature,  T^  is  given  by  the  equation: 
Tj       =      T.  +  PdGj,  (eqn3.7) 
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Figure  3.3    P-I-N  Diode  C  vs.  I-Region  Diameter  and  Thickness. 
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Figure  3.4    Junction  Voltage  at  High  RF  and  Reverse  Bias. 
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The  value  of  0:^  is  determined  by  the  physical  volume  of  the  junction, 
material  and  geometry  of  the  package.    The  maximum  value  of   9:^  for  most  glass 
packages  can  be  as  much  as  600°C/watt.   A  large-area,  high  capacitance  chip  mounted 
on  copper  heat  sink  can  have  extremely  low    0-    of  5°C/watt  or  less  [Ref  9:  p.  10]. 
b.    Pulse  Power  Dissipation 

For  a  short  duration  pulse,  the  temperature  limit  is  determined  by  the 
thermal  time  constant  x^-^  o^  the  active  region. 

T^^      =      e-HC  (eqn3.8) 

where    0:  is  the  junction  region  thermal  resistance,  this  value  is  available  in  the 
specification  sheet  or  can  be  calculated  from  the  diode  physical  size. 


W 
0.        =      — -  (eqn3.9) 

where  W     =    the  "width  of  the  I-region, 

K^j^  =    the  thermal  conductivity  of  the  material, 

A      =    the  junction  area,  and 

HC    =     the   thermal   capacity    of  the   active   region,    also    available   in   the 
specification  sheet  or  can  be  calculated  from: 

HC       =      C  pWA  (eqn3.10) 

where  C_    =    the  specific  heat  of  the  material,  and 
p      =    the  density  of  the  material. 

For  a  pulse  width  of  0.25  t^j^  or  less,  the  pulse  width  is  small  compared 
with  the  thermal  time  constant.   All  the  thermal  energy  is  absorbed  in  the  active  region 

area,  not  spreading    to  the  surrounding.  The  temperature  nse  of  the  junction.  AT,  m 
units  of       C  is  exoressed  bv: 


Ph  At 
AT       =     -S —  (eqn3.ll) 

HC 
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where  At  is  the  pulse  duration.  The  junction  temperature^  T:  for  a  short  pulse  width 
from  the  above  assumption  is  the  sum  of  the  ambient  temperature  and  the  temperature 
rise  of  the  junction. 


T-       =     T^  +  AT  (eqn3.12) 

For  pulse  widths  of  0.50t^j^  to  3t^j^  the  diode  junction  temperature  follows 
an  exponential  curve  given  by  the  equation 

T|       =     T^  +  Pd  ej  [1  -  exp(-t/T^i^)  ]  (eqn  3.13) 

Note  that  equation  3.13  contain  0-,  the  junction  thermal  resistance  and  not  0:„  the 
total  thermal  resistance.  For  cooling  between  pulses,  the  junction  temperature  follows 
an  exponential  decay  with  the  same  time  constant  t^-^  [Ref.  9:  p.  10]. 

D.      TABLE  OF  TYPICALLY  AVAILABLE  P-I-N  DIODES 

To  illustrate  typical  p-i-n  diodes  available  for  microwave  switch  design,  Table  1 
lists  a  number  of  typical  diode  chips  from  Microwave  Associates,  Inc.  and  their 
approximate  electrical  parameters.  The  symbols  used  in  this  table  are  already  defined  in 
this  chapter. 


^Neglecting  the  heat  fiow  from  the  junction. 
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TABLE  1 

TYPICAL  PARAMETERS  OF  P-I-N  DIODES 

MODEL 

MA-        MA- 

MA- 

MA- 

MA-        MA- 

MA- 

NUMBER 

4P102     4P103 

4P202 

4P203 

4P303      4P404 

4P504 

Minimum 

Vb(VoIt) 

50             30 

100 

100 

200          300 

500 

Maximum 

Cj  @  Vr(lO) 

0.5            .15 

.05 

.15 

.15            .20 

.20 

(pF)(Volt) 

Maximum 

Rs  @  IfllO) 

2.5            1.5 

2.0 

1.5 

1.5           0.6 

0.6 

(ohm){ma) 

Typical 

Carrier 

.100          .010 

.100 

.100 

.200           1.00 

2.00 

Liretime(jLiS) 

Maximum 

Thermal 

60             40 

60 

30 

30          20 

20 

Resistance 

("C/watt) 

Power  (max) 

Dissipation 

2.5            3.75 

2.5 

5.0 

5.0         7.5 

7.5 

at25°C 

(watt) 
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IV.  EVALUATION  OF  SEMICONDUCTOR  MICROWAVE  SWITCH 

CHARACTERISTICS 

In  this  chapter,  semiconductor  device  figures  of  merit,  Kawakami's  M  and  diode 
Q,  are  tested  for  the  shunted  switch  configuration,  to  find  the  best  parameter  for  use  in 
a  s"witch  selection  and/or  design  problem.  From  equation  1.5,  we  know  that  the  figure 
of  merit  M,  is  a  function  of  insertion  loss,  S2|(l),  in  the  switch  on-state,  and  load 
isolation,  S2j{2),  in  the  switch  off-state.  In  order  to  solve  the  equation  mentioned 
above,  the  switching  characteristics,  insertion  loss  and  load  isolation,  should  be 
provided.  By  using  a  CAD  program  (Touchstone)  we  can  minimize  the  insertion  loss  in 
the  switch  on-state,  and  maximize  the  load  isolation  in  its  off-state.  From  the  value  of 
these  two  switching  characteristics,  we  can  compare  and  select  the  figure  of  merit  that 
best  characterizes  the  performance  of  this  type  of  microwave  switch.  The  power 
dissipation  and  heat  removal  requirements  in  the  individual  semiconductor  device  can 
be  estimated  from  its  terminal  voltage.  We  will  introduce  two  methods  to  evaluate  the 
semiconductor  device  terminal  voltage,  and  its  power  dissipation: 

1.  Using  a  Touchstone  application  program  to  evaluate  the  terminal  voltage  of 
interest. 

2,  Using  a  Fortran  program  written  from  the  circuit  model  and  method  introduced 
in  Chapter  II. 

A.       INSERTION    LOSS    AND    LOAD    ISOLATION    CALCULATION    USING 
COMPUTER-AIDED  CIRCUIT  ANALYSIS  PROGRAM  (TOUCHSTONE) 

The  Touchstone  program  is  advanced  software  for  RF/microwave  computer- 
aided  engineering  [Ref  11].  It  can  provide  an  output  of  the  S-parameters  of  a  multi- 
port  network,  and  can  perform  an  interactive  optimization,  and  has  applications  which 
can  be  used  in  the  measurement  of  selected  terminal  voltages.  For  a  given  electrical 
circuit,  we  can  create  a  circuit  file  and  evaluate  circuit  performance  or  use  the 
optimization  mode  to  get  the  optimized  circuit-component  values. 

There  are  two  alternative  methods  of  approach  to  analyze  these  two  switching 
characteristics.  The  first  method  uses  semiconductor  device  parameters,  diode  junction 
resistance  (R:)  and  diode  junction  capacitance  (C;),  as  the  direct  input  parameters  of 
the  circuit  file.  The  second  method  employs  semiconductor  device  figure  of  merit,  diode 
Q,  as  an  input  parameter  of  the  circuit  file.    In  this  approach  the  diode  junction 
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capacitance  and/ or  resistance  are  calculated  in  the  circuit  file,  and  then  follow  the  same 
procedure  as  the  first  category. 

Two  types  of  microwave  switches,  single  pole  single  throw,  and  single  pole 
double  throw  switches,  with  one  or  multiple-semiconductor  devices  are  investigated.  By 
using  the  electrical  equivalent  circuit  model  introduced  in  Chapter  II,  we  can  create  the 
circuit  file  following  one  of  the  above  categories  for  the  Touchstone  system.  The 
assumptions  for  the  circuit  model  are: 

1.  All  diodes  in  the  switch  are  identical. 

2.  The  transmission  line  configurations  in  the  switch  are  symmetrical. ' 

3.  The  transmission  line  sections  are  lossless. 

4.  The  transmission  line  impedance  is  between  30  to  90  Ohms. 

5.  The  transmission  line  electrical  length  is  between  30  to  120  degrees. 

If  the  input  of  the  circuit  file  is  the  figure  of  merit,  or  diode  Q,  then  the  junction 
resistance  and  junction  capacitance  are  calculated  from  equations  4.2  and  4.3: 


1000 


1000 
R-  =         {eqn4.2) 

2^FoCiQdiode 


J 


1000 
C;  =  (eqn4.3) 

2^Po<^iQdiode 


'J 


where  Q^i^ode  ^^  ^^^  figure  of  merit. 

R:  is  the  diode  junction  resistance  in  ohm. 
C;  is  the  diode  junction  capacitance  in  pF. 
Fq  is  the  center  frequency  in  GHz. 

The  relation  of  R-,  C-  and  diode  Q  at  the  center  frequency  of  15  GHz.  from  equation 
4.2,  4.3  is  shown  m  Figure  4.1.  In  our  circuit  file  for  convenience  we  fix  the  value  of 
capacitance  and  calculate  the  value  of  junction  resistance  (equation  4.2)  and  var\'  the 
value  of  diode  Q. 


^For  two-diode  switches,  TLINl  is  identical  to  TLIN3  (see  Figures  4.3  and  4.7). 
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Figure  4.1     Resistance  vs  Capacitance  as  a  Function  of  Diode  Q. 
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In  the  procedure  of  creating  a  circuit  file  to  optimize  the  switch  characteristics 
for  low  insertion  loss  and  high  isolation,  one  more  assumption  is  made;  the  impedance 
values  of  the  diode  are  kept  unchanged  at  the  particular  frequency  in  each  state.  The 
length  and  the  characteristic  impedance  of  the  transmission  line  sections  are  varied  to 
achieve  the  optimum  switching  characteristics  in  the  operating  frequency  interval.  In 
the  optimization  mode  of  the  circuit  file,  we  can  minimize  the  switch  insertion  loss  and 
maximize  load  isolation  at  the  same  time.  This  approach  may  require  a  long 
computation  time  to  produce  a  result,  and  the  switch  insertion  loss  still  may  not 
approach  the  wanted  value.  It  is  found  that  a  more  efficient  method  of  approach  is  to 
optimize  only  the  insertion  loss.  The  computation  time  is  typically  shorter,  and 
normally  produces  satisfactory  insertion  loss  and  load  isolation.  Figure  4.2  is  the  flow 
chart  of  steps  in  creating  a  circuit  file  for  the  Touchstone  system. 

1.  Insertion  Loss  in  the  Snitch  On-state 

In  the  switch  on-state,  the  ratio  of  the  available  power  to  the  power  reaching 
the  load  is  defined  as  the  insertion  loss  of  the  switch.  In  this  switching  state,  the  p-i-n 
diode  is  equivalent  to  a  very  high  impedance  in  shunt  with  the  circuit.  A  small  amount 
of  energy  will  be  dissipated  in  the  diode,  but  most  of  the  energy  is  delivered  to  the  load. 
Therefore  the  insertion  loss  of  a  well  designed  switch  should  be  as  low  as  possible.  In 
order  to  achieve  the  low  insertion  loss,  the  optimization  mode  in  the  Touchstone 
system  is  used  to  perform  the  minimization  of  insertion  loss  to  obtain  the  circuit  design 
values.  In  practice  the  insertion  loss  has  units  of  dB,  and  can  be  calculated  from: 

[p.  ] 
Insertion  loss  in  dB     =         lOlog^Q -^^^^ —  (eqn  4.4) 

t^load 

From  Chapter  II,  under  the  assumption  of  an  ideal  switch  with  matched  load  at  both 
input  and  output  ports,  we  can  use  the  scattering  matrices  to  represent  the 
transmission  behavior  for  both  types  of  switches.  Recall  from  Chapter  II,  for  the  on- 
state  reDresentadon  of  S-matrices: 


SPST  ON  STATE:   [S] 


0  1 

1  0 


(eqn  4.5) 
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SPDT  ARM  (2)  ON:     [S]    = 


0 

1 

0 

1 

0 

0 

0 

0 

0 

(eqn  4.6) 


From  equations  4.5  and  4.6,  the  insertion  loss  of  the  switch  is  defined  as  !S-.|I  in  dB^ 
and  can  be  evaluated  from  the  circuit  file  of  Touchstone  system. 
2.  Load  Isolation  in  the  S>vitch  Off-state 

In  the  switch  off-state,  the  load  isolation  of  the  switch  is  also  defined  by  the 
ratio  of  the  available  power  to  the  power  reaching  the  load.  In  this  switching  state,  the 
p-i-n  diode  is  equivalent  to  a  very  low  impedance  in  shunt  with  the  circuit.  Most  of  the 
energy  is  reflected  back  to  the  source,  some  of  the  energy  is  dissipated  in  the  diode  and 
a  small  amount  of  energy  can  transfer  to  the  load.  For  an  ideal  switch  in  its  off-state, 
the  load  isolation  should  be  infinite,  meaning  that  no  energy  is  transferred  to  the  load. 
In  the  real  switch,  the  higher  load  isolation  is  the  better  switch  performance.  The  load 
isolation  has  units  of  dB  and  has  the  same  formula  as  in  equation  4.4.  In  order  to 
achieve  the  best  performance  of  the  designed  swiT:ch,  the  Touchstone  optimization 
mode  is  used,  if  desired,  to  optimize  the  load  isolation  as  well  as  insertion  loss.  The  S- 
matrix  representation  of  the  switch  in  the  off-state  with  the  same  assumptions  as  for 
the  insertion  loss  is  defmed  as  follows: 


SPST  OFF  STATE:  [S]      = 


1    0 
0     1 


(eqn  4.7) 


SPDT  ARM  (3)  ON:    [S]    = 


0 

0 

1 ' 

0 

0 

0 

1 

0 

0 

(eqn  4.S) 


From  [Refs.  5.61,  and  equations  4.7  and  4.8  IS^^il  in  dB  is  denned  as  the  load  isolation 
of  the  SPST  switches,  and  jS^^i  m  dB  is  defined  as  the  load  isolation  oi  the  SPDT 
switches. 


^dB  value  =  20  log^Q  S21 
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INPUT  PARAMETERS: 

1.  TRANSMISSION  LINE 
Initial  Characteristic 
Impedance 

Initial  Zlecrrical  Length 

2.  DIODE  PARAMETERS 

Junction  Capacitance, Resistance 

3.  OPERATING  FREQUENCY 


V/L 


OUTPUTS: 
SWITCH  CHARACTERISTICS 
Inser1?ion  LOss 
Load  Isolation 


MINIMIZE   INSERTION  LOSS 

1.  TRANSMISSION  LINE 
Characteristic  Impedance 
Electrical  Length 

2.  SWITCH  CHARACTERISTICS 
Insertion  Loss 

Load  Isolation 


io 


-^TOPj 


Figure  4.2     Flow  Chart  of  a  Circuit  File. 
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3.  Example  of  Insertion  Loss  and  Load  Isolation  Evaluation 

We  divide  the  evaluation  of  switciiing  characteristics  into  two  cases: 

1.  Single  pole  single  throw  switch  with  two  diodes. 

2.  Single  pole  double  throw  switch  with  four  diodes.^ 

Two  categories  of  creating  the  circuit  file  are  applied  to  both  cases. 
a.   SPST  Two- Diode  Switch 

The  electrical  circuit  representation  of  the  SPST  two-diode  switch  is  shown 
in  Figure  4.3(a).  The  switch  on-  and  off-state  equivalent  circuits  are  shown  in  Figures 
4.3  (b)  and  (c)  respectively.  To  create  a  circuit  file  from  the  equivalent  circuit,  node 
numbers  are  assigned  to  every  component,  then  the  steps  in  Figure  4.2  are  followed. 
The  example  of  circuit  files  for  switching  characteristics,  insertion  loss  and  load 
isolation,  evaluation  of  a  SPST  two-diode  switch  are  listed  in  Appendix  A. 

SPSTD.CKT,  listed  in  Appendix  Al,  is  an  example  of  the  circuit  file 
created  for  the  single  pole  single  throw  switch  in  the  first  category,  in  which  diode 
resistance  and  capacitance  values  are  chosen  and  used  as  input  data  for  the  circuit  file. 
Input  parameters  are  diode  junction  resistance  and  diode  junction  capacitance  from 
typical  p-i-n  diode  parameters  available  in  Table  1.  We  select  the  diode  model  number 
MA-4P203  with  0.15  pF  junction  capacitance  and  1.5  ohms  junction  resistance  to  be 
the  input  parameters  for  the  circuit  file.  After  optimizing  the  circuit  file  from  10  to  20 
GHz,  we  have  the  output  file  of  the  switch  insertion  loss  and  load  isolation  as  shown 
in  Table  2. 

From  the  Touchstone  output  file  listing  in  Table  2,  LINDON  and 
LINDOFF  are  switch  insertion  loss  and  load  isolation  respectively.  Our  required 
switch  operating  frequency  range  is  between  10  and  20  GHz,  the  wanted  value  of 
switch  insertion  loss  is  around  -0.1  to  -0.2  dB  and  the  load  isolation  greater  than  -45 
dB.  The  characteristics  of  Table  2  satisfy  our  switch  design  requirement.  The  plot  of 
insertion  loss  and  load  isolation  is  shown  in  Figure  4.4. 

After  optimization,  the  circuit-component  values,  characteristic  impedances 
and  electrical  lengths  of  the  transmission  line  sections,  are  given  in  the  optimized 
circuit  file.  The  optimized  line-section  parameters  of  the  example  m  Appendix  Ai  are: 
transmission  line  section  A  has  Z^  of  39.40  ohms  and  E^  of  106.3  degrees,  transmission 
line  section  B  has  Z^  of  38.25  ohms  and  E^  of  62.7  degrees.  In  switch  construction,  to 
obtain  the  design  value  of  the  insertion  loss  and  load  isolation,  the  transmission  line 


^There  are  two  diodes  in  each  branch  of  the  switch. 
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4.3(a)  SPST  2-DIODE  EQUIVALENT  CIRCUIT 
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4.3(b)  SPST  ON-STATE  EQUIVALENT  CIRCUIT 


TLIN  A   2    TLIN  B    3    TLIN  A 
1  1      ^ 1      I ^ 1      I 4 


Rj 


Rj 


4.3(C)  SPST  OFF-STATE  EQUIVALENT  CIRCUI' 


Figure  4.3     Single  Pole  Single  Throw  Switch.  Equivalent  Circuit. 
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Figure  4.4    SPST  Switch  Insenion  Loss  and  Isolation  vs  Frequency. 
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TABLE  2 

TOUCHSTONE  OUTPUT  SPSTD.CKT 

Frequency  - 

GHz 

dB  [S21] 

dB  [S21] 

LINDON 

LINDOFF 

10.0000 

-0.150 

-46.276 

12.0000 

-0.053 

-47.025 

14.0000 

-0.085 

-48.011 

16.0000 

-0.110 

-49.261 

18.0000 

-0.141 

-50.590 

20.0000 

-0.195 

-51.750 

sections  in  the  switch  should  have  their  characteristic  impedances  and  electrical  lengths 
close  to  the  calculated  values. 

SPSTQ.CKT,  listed  in  Appendix  A2,  is  an  example  of  a  circuit  file  created 
for  the  single  pole  single  throw  switch  in  which  input  parameters  are  Q^j^ode  ^^^  diode 
junction  capacitance.  The  selected  value  of  junction  capacitance  is  0.15  pF  and  the 
value  of  diode  Q  are  between  10  and  500  are  used  as  input  parameters  of  this  circuit 
file.  After  optimizing  the  circuit  file  from  10  to  20  GHz  we  have  the  output  file  of  the 
switch  insertion  loss  and  load  isolation  as  shown  in  Table  3. 

Table  3  is  the  collection  of  of  SPSTQ.CKT  outputs  evaluated  for  different 
diode  Q  values.  This  table  verifies  that  the  diode  Q  has  the  character  of  a  figure  of 
merit,  serving  as  an  indicator  of  switching  performance.  The  higher  diode  Q  is,  the 
better  the  suitcii  characteristics:  low  insertion  loss  and  high  isolation. 

Figures  4.5  and  4.6  show  cne  plot  of  insertion  loss  vs  frequency  and  load 
isolation  vs  frequency  of  a  SPST  switch  with  two  diodes  with  diode  Q  as  a  parameter. 
We  can  use  these  two  figures  together  with  Figure  4.1  for  a  design  problem  of  a  SPST 
switch  from  a  given  specification.  For  example,  a  specification  may  require  a  two  p-i-n 
diode  SPST  switch  to  have  at  least  -0.1  dB  insertion  loss  and  isolation  better  than 
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TABLE  3 

TOUCHSTONE  OUTPUT  SPSTQ.CKT 

INSERTION  LOSS  IN  dB 

Frequency 

Diode 

Diode 

Diode 

Diode 

Diode 

Diode 

GHz 

Q=10 

Q=20 

Q=50 

Q=100 

Q  =  200 

Q=500 

10.0000 

-.367 

-.342 

* 

-.135 

-.064 

-.056 

-.052 

12.0000 

-.194 

-.111 

-.049 

-.031 

-.019 

-.011 

14.0000 

-.342 

-.170 

-.082 

-.062 

-.044 

-.034 

16.0000 

-.499 

-.231 

-.101 

-.052 

-.037 

-.013 

18.0000 

-.517 

-.269 

-.134 

-.072 

-.039 

-.019 

20.0000 

-.724 

-.393 

-.186 

-.096 

-.048 

-.019 

LOAD  ISOLATION  IN  dB 

Frequency 

Diode 

Diode 

Diode 

Diode 

Diode 

Diode 

GHz 

Q=10 

Q=20 

Q  =  50 

Q=100 

Q=200 

Q  =  500 

10.0000 

-20.21 

-30.53 

-47.38 

-58.78 

-70.62 

-86.42 

12.0000 

-20.51 

-30.84 

-48.06 

-59.46 

-71.30 

-87.09 

14.0000 

-21.05 

-31.37 

-48.88 

-60.53 

-72.37 

-88.17 

16.0000 

-21.89 

-32.27 

-50.17 

-61.90 

-73.76 

-89.57 

18.0000 

-22.93 

-33.42 

-51.39 

-63.26 

-75.14 

-90.96 

20.0000 

i 

-23.95 

-34.57 

-52.40 

-64.28 

-76.13 

-92.01 

-60.00  dB  at  the  mid- frequency  of  15  GHz.  In  Figure  4.5,  with  the  value  of  diode  Q 
50  or  more  the  switch  can  have  insertion  loss  less  than  -0.1  dB.  But  from  Figure  4.6 
the  diode  Q  of  50  is  not  great  enough  to   obtain  load  isolation  of -60.00  dB. 
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Figure  4.5     SPST  Switch  Insertion  Loss  vs  Q^^jode* 
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Figure  4.6    SPST  Switch  Load  Isolation  vs  Q^jjode" 


47 


Therefore  we  select  a  diode  for  the  switch  that  has  the  value  of  Q<;iiQ(je  of  100. 
From  Figure  4.1  with  the  Q=  100  curve  we  can  choose  one  of  the  p-i-n  diodes  that 
has  the  value  of  junction  capacitance  and  resistance  close  to  our  estimated  value.  If 
we  select  p-i-n  diode  model  number  MA-4P504  with  diode  junction  capacitance  of  0.2 
pF  and  junction  resistance  of  0.6  ohm,  at  15  GHz  we  obtain  the  diode  Q  value  o'C 
88.42,  switch  insertion  loss  of  -.102  dB  and  isolation  of  -60  dB  from  the  modified 
SPST.CKT  circuit  file. 

b.  SPOT  Four-Diode  Switch 

A  single  pole  double  throw  switch  is  considered  as  two  SPST.switches  with 
input  arms  tied  together,  when  one  arm  is  in  the  on-state,  the  other  is  in  the  off-state. 
A  SPDT  four-diode  switch  has  the  electrical  circuit  representation  and  the  switch  on 
off-state  equivalent  circuit  shown  in  Figure  4.7(a)  and  (b)  respectively.  To  create  a 
circuit  file  from  the  equivalent  circuit,  the  same  procedure  can  be  followed  as  in  SPST 
switch  case.  Or  we  can  make  some  necessary  modifications  to  the  SPST  switch  circuit 
files  in  Appendix  A.  The  example  of  circuit  files  for  evaluating  SPDT  switch  insertion 
loss  and  load  isolation  are  listed  in  Appendix  B. 

SPDTD.CKT,  listed  in  Appendix  Bl,  is  an  example  of  the  circuit  file  used 
for  the  single  pole  double  throw  switch  with  four  diodes.  With  the  same  p-i-n  diode 
model  number  as  in  the  circuit  file  SPSTD.CKT,  after  optimizing  this  circuit  file  from 
10  to  20  GHz,  we  have  the  output  file  of  the  insertion  loss  and  load  isolation  shown  in 
Table  4.  The  optimized  transmission  line  section  characteristic  impedances  and 
electrical  lengths  are:  46.56  ohms  and  93.58  degrees  for  transmission  line  section  A, 
89.06  ohms  and  32.82  degrees  for  transmission  line  section  B. 

The  output  file  listing  in  Table  4  shows  that  for  the  required  operating 
frequency  between  10  to  20  GHz,  the  designed  SPDT  four-diode  switch  using  p-i-n 
diode  model  number  MA-4P203  with  0.15  pF  junction  capacitance  and  1.5  ohms 
junction  resistance,  has  insertion  loss  around  -1.0  dB  and  load  isolation  greater  than 
-50  dB  We  realize  that  the  insertion  loss  of  the  SPDT  switch  case  is  not  so  low  as  for 
the  SPST  switch  case.  However  the  insertion  loss  and  load  isolation  satisty  the  general 
requirement  of  a  switch  design  problem.  The  plot  of  insertion  loss  and  load  isolation 
vs  operating  frequency  is  shown  in  Figure  4.8. 

SPDTQ.CKT,  listed  in  Appendix  B2,  is  an  example  of  the  circuit  file 
created  for  the  single  pole  double  throw  switch  for  which  input  parameters  are  Q(jiode 
and  diode  junction  capacitance.     In  this  example,   the  value   of  0.15  pF  junction 
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TABLE  4 

TOUCHSTONE  OUTPUT  SPDTD.CKT 

Frequency  - 

GHz 

dB  [S21] 

dB  [S21] 

LINDON 

LINDOFF 

10.0000 

-0.756 

-56.199 

12.0000 

-0.327 

-57.856 

14.0000 

-0.226 

-59.898 

16.0000 

-0.385 

-61.603 

18.0000 

-0.417 

-61.471 

20.0000 

-0.665 

-59.573 

capacitance,  and  values  of  diode  Q  from  10  to  500  are  selected  and  used  as  input 
parameters  of  the  circuit  file.  The  optimized  output  file  SPDTQ.CKT  is  shown  in 
Table  5. 

For  general  use  of  diode  Q  as  the  figure  of  merit  to  characterize  the 
switching  performance  of  this  type  of  switch,  the  choice  of  junction  capacitance  value 
selected  as  input  parameter  to  the  Touchstone  file  does  not  strongly  affect  the  Q- 
dependence  of  switch  performance  shown  in  Table  5.  The  result  of  output  file 
SPDTQ.CKT  listed  in  Table  5  verifies  that  the  diode  Q  has  the  character  of  a  figure  of 
merit,  serving  as  an  indicator  of  the  switch  performance.  In  a  switch  design  problem, 
diode  Q  can  be  used  as  a  parameter  that  characterizes  the  switching  transmission 
characteristics.  The  plot  of  insertion  loss  and  load  isolation  vs  operating  frequencv  as 
a  function  of  diode  Q  from  the  results  m  Table  5  are  shown  m  Figures  4.9  and  4.10. 
These  two  plots  can  be  used  m  a  design  problem  of  a  SPDT  switch  iroin  a  given 
specification.  The  procedures  for  the  design  of  a  SPDT  switch  are  the  same  as 
mentioned  in  the  SPST  switch  case. 
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TABLE  5 

TOUCHSTONE  OUTPUT  SPDTQ.CKT 

INSERTION  LOSS 

INdB 

Frequency 

Diode 

Diode 

Diode 

Diode 

Diode 

Diode 

GHz 

Q=10 

Q  =  20 

Q  =  50 

Q=100 

Q  =  200 

Q  =  500 

10.0000 

-1.96 

-1.38 

-1.04 

-.926 

-.824 

-.796 

12.0000 

-1.36 

-.725 

-.345 

-.217 

-.148 

-.117 

14.0000 

-1.13 

-.577 

-.265 

-.179 

-.129 

-.095 

16.0000 

-1.32 

-.756 

-.394 

-.294 

-.232 

-.190 

18.0000 

-1.69 

-.935 

-.395 

-.214 

-.118 

-.059 

20.0000 

-2.37 

-1.49 

-.835 

-.564 

-.515 

-.425 

LOAD  ISOLATION  IN  dB 

Frequency 

Diode 

Diode 

Diode 

Diode 

Diode 

Diode 

GHz 

Q=10 

Q  =  20 

Q=50 

Q=100 

Q  =  200 

Q  =  500 

10.0000 

-33.55 

-44.74 

-59.50 

-71.18 

-83.59 

-99.46 

12.0000 

-34.35 

-45.92 

-60.83 

-72.58 

-85.09 

-100.95 

14.0000 

-35.82 

-47.90 

-63.08 

-74.99 

-87.56 

-103.43 

16.0000 

-37.56 

-49.71 

-64.75 

-76.66 

-89.13 

-105.00 

18.0000 

-38.27 

-49.68 

-64.27 

-75.97 

-88.18 

-105.04 

20.0000 

-37.54 

-48.24 

-62.75 

-74.37 

-86.21 

-102.20 

B.       THE   SEMICONDUCTOR   DEVICE    FIGURE    OF   MERIT   WHICH   BEST 
CHARACTERIZES  THE  SWITCHING  PERFORMANCE 

Two  well  known  figures  of  merit  were  used  as  parameters  to  predict  the  switching 

performance    of  semiconductor   microwave    switches.     In   this    thesis    research,    the 
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Figure  4.10    SPDT  Switch  Load  Isolation  vs  Q^Q^jg. 
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semiconductor  device  considered  is  tlie  p-i-n  diode,  used  in  a  sliunt-configuration 
microwave  switch.  The  conventional  diode  Q  was  shown  to  coincide  with  the 
simplified  equation  of  Kurokawa  and  Schlosser  for  quality  factor  Q,  and  therefore  we 
use  the  conventional  Q  as  our  semiconductor  device  figure  of  merit.  In  Chapter  II 
equations  2.10  and  2.11  showed  the  relationship  between  both  figures  of  merit:  the 
diode  Q  and  Kawakami's  M.  To  verify  these  two  equations,  we  investigated  the 
relationship  of  Kawakami's  M  and  the  switch  transmission  characteristics  defined  in 
equation  1.5.  Optimal  values  of  switch  insertion  loss  and  load  isolation  were  evaluated 
from  the  Touchstone  circuit  file  which  has  the  diode  Q  as  an  input  parameter.  For  the 
SPDT  s'W'itch  case,  values  of  insertion  loss  and  load  isolation  evaluated  with  different 
diode  Q  values  at  15  GHz.  are  substituted  in  equation  1.5,  and  the  results  are  listed  in 
Table  6. 


TABLE  6 

COMPARISON  BETWEEN  FIGURES  OF  MERIT  M, 

AND  DIODE  Q 

Diode  Q 

S21  (1) 

S21  (2) 

M  (eqn  1.5) 

AM 

—  xi00(%) 

AQ 

10.0000 

0.9421413 

0.1590319 

0.9770399 

0.1667 

20.0000 

0.9677804 

0.0865225 

0.9937131 

0.0161 

50.0000 

0.9830945 

0.0404859 

0.9985498 

1.8e-3 

100.000 

0.9872322 

0.0222759 

0.9994510 

2.9e-4 

200.000 

0.9897033 

0.0118343 

0.9997507 

5.2e-5 

500.000 

0.9915362 

0.0053798 

0.9999076 

Table  6    shows  the   comparison    of  the  two  figures  of  merit,  with  M    calculated 
from    equation  1.5.    Note    that  in    this    table    S2i(l),  S2j(2)    correspond    to  switch 
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insertion  loss  and  load  isolation  respectively,  where  the  value  of  S2J  is  in  the 
form:  $21    =    10^*^°'-^^'.    The  values  of  M  are  less  than  1  as  mentioned  earlier.    For 

this  type  of  switch  at  Q^jjode  °^  ^^  ^^^  "^^^^^  °^  ^^  ^^  0.9762,  and  at  Q^j^ode  °^500  the 
value  of  M  is  0.9999.  For  typical  diode  Q  values,  M  departs  very  little  from  unity,  and 
the  percentage  of  deviation  is  reduced  tremendously  in  the  high  Q  region.  Thus  the 
parameter  Q^jjode  ^^  niore  convenient  to  use  due  to  its  wider  range  compared  with  the 
parameter  M.  When  tested  for  different  types  of  switch, ^^  the  results  are  nearly  the 
same  as  shown  in  Table  6  for  all  the  values  of  M  (almost  unity  and  small  deviation  in 
each  interval).  An  additional  reason  in  support  of  our  decision  to  select  the 
conventional  diode  Q  as  the  best  figure  of  merit  to  characterize  switching  performance 
is  that  it  is  well  kno\^Ti  and  widely  used.  Table  7  shows  the  comparison  of  the  result 
from  Table  6  and  values  of  M  calculated  from  equations  2.10  and  2.11,  this  table 
verifies  that  equation  2.10  can  be  used  to  convert  the  diode  Q  value  to  the  figure  of 
merit  iVI,  and  vice  versa. 


-■   ■ 

TABLE  7 

COMPARISON  OF  FIGURES  OF  MERIT  M, 

APPROACHES 

FROM  DIFFERENT 

Diode  Q 

M  (eqn  1.5) 

M  (eqn  2.10) 

M  (eqn  2.11) 

10.0000 

0.9770399 

0.9805806 

0.9800000 

20.0000 

0.9937131 

0.9950372 

0.9950000 

50.0000 

0.9985498 

0.9992001 

0.9992000 

100.0000 

0.9994108 

0.9998000 

0.9998000 

200.0000 

0.9997507 

0.9999499 

0.9999500 

500.0000 

0.9999076 

0.9999919 

0.9999920 

1 

10 


SPDT  s\\itch  with  tuning  stub,  micro  strip  model. 
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C.       NODAL  VOLTAGE  CALCULATION  USING  TOUCHSTONE  PROGRAM 

The  power  dissipation  in  a  particular  diode  can  be  calculated  by  the  voltage 
across  its  resistance  or  the  current  flow  in  the  resistor  of  the  p-i-n  diode  equivalent 
circuit.  In  the  Touchstone  program  there  is  a  technique  for  probing  the  nodes  within  a 
circuit  file  to  determine  the  voltages  at  those  nodes.  An  ideal  voltage-controlled 
voltage  source  (VCVS)  is  used  as  a  non-loading  voltmeter  probe  [Ref  12].  The  input 
terminals  of  the  VCVS  are  connected  to  the  nodes  between  which  the  voltage  is  to  be 
determined.  The  output  terminals  are  connected  to  an  output  port  so  that  the  voltages 
may  be  read.  The  voltages  read  out  from  this  method  are  normalized  to  unity  value  of 
the  open  circuit  generator  voltage  [Ref  12].  Appendix  C  shows  examples  of  the 
Touchstone  circuit  files  written  for  the  evaluation  of  nodal  voltages  of  SPST  and  SPDT 
switches.  For  SPDT  switch  case,  the  equivalent  circuit  model  using  in  nodal  voltage 
calculation  is  shown  in  Figure  4.11. 

The  circuit  file  SPDTV.CKT,  listed  in  Appendix  C2,  is  an  example  of  a 
Touchstone  circuit  file  using  VCVS  modules  acting  as  voltage  probes  reading  the  nodal 
voltages  from  the  attached  terminals.  The  input  parameters  for  this  circuit  file  are  the 
diode  parameters  and  the  optimal  transmission  line  section  values  of  the  SPDT  switch 
evaluated  in  Appendix  B2.  The  major  disadvantage  of  the  Touchstone  circuit  file  is 
that  we  can  define  up  to  four  ports  in  our  circuit  file,  meaning  that  we  can  measure  no 
more  than  three  nodes  at  a  time.  To  measure  a  SPDT  switch  with  four  diodes,  there 
are  six  nodes  of  interest.  The  first  three  nodes  are  measured  in  the  SPDTV.CKT  circuit 
file,  while  the  last  three  nodes  are  measured  by  making  some  modifications  to  this 
circuit  file  (by  moving  the  VCVS  module  to  the  new  positions).  The  output  files  of 
SPDTV.CKT  are  listed  in  Table  8. 

Knowing  terminal  voltage  across  the  resistive  part  of  the  p-i-n  diode  equivalent 
circuit  both  on  and  off-state,  the  next  step  is  to  calculate  the  power  dissipation  in  the 
individual  p-i-n  diode  from  the  equation: 

P^     =         -^  (eqn4.9) 

where    P^^  is  the  dissipated  power  at  node  n, 

y^  is  the  terminal  voltage  at  node  n,  and 

Rj^  is  the  resistive  component  of  the  n^^  p-i-n  diode. 
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Figure  4.11     Nodal  Voltage  Measurement  Equivalent  Circuit. 


TABLE  8 

NODAL  VOLTAGE  FROM  TOUCHSTONE  SPDTV.CKT 

SWITCH  ON-STATE 

SWITCH  OFF-STATE 

Frequency 

^load 

Vdi 

Vd2 

^load 

Vdi 

Vd2 

10.0000 

.458 

.006 

.006 

7.7e-4 

.017 

7.3e-4 

12.0000 

.482 

.008 

.008 

6.4e-4 

.016 

6.0e-4 

14.0000 

.485 

.010 

.009 

5.1e-4 

.015 

4.7e-4 

16.0000 

.478 

.012 

.010 

-  2e— i 

.014 

3.9e-^ 

18.0000 

.477 

.015 

.011 

4.2e-4 

.015 

4.0e-4 

20.0000 

.463 

.018 

.012 

5.3e-4 

.021 

5.0e-4 
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In  practice  the  dissipated  power  in  a  particular  diode  is  normalized  to  the  power 
delivered  to  the  load  in  the  switch  on-state.  Therefore  each  power  dissipation  calculated 
from  equation  4.9  is  normalized  to  the  switch  on-state  output  power.  For  the 
normalized  power  value,  equation  4.9  becomes: 

p^    =  ^^n^    ^joad  (eqn4.10)    * 

f^loaJ    ^n 

where  V^^^^  is  the  load  terminal  voltage  of  the  switch  on-state. 

The  use  of  Touchstone  circuit  files,  SPSTV.CKT  and  SPDTV.CKT,  is  very 
convenient  for  evaluating  the  diode  terminal  voltage,  but  we  cannot  use  Touchstone  to 
solve  for  power  dissipation  and  the  diode  junction  temperature  problem  directly.  To  do 
so  one  may  use  a  hand  calculator  or  ""vrite  a  special  computer  program  calculating  the 
power  dissipation  in  each  diode  from  equation  4.10,  and  use  these  dissipated  powers  to 
calculate  the  heat  production  from  equations  in  Chapter  III.  However  we  see  that  the 
Touchstone  output  file  has  only  three  significant  digits  after  the  decimal  point, 
meaning  that  numerical  round  off  takes  place  in  the  third  digit.  These  outputs 
generally  satisfy  engineering  requirements.  However  for  a  very  small  nodal  voltage 
across  the  diode,  these  rounded  off  outputs  may  cause  significant  error  when  we 
substitute  them  in  a  quadratic  equation  such  as  equation  4.10.  It  is  a  serious 
disadvantage  to  use  these  outputs  for  high  output  power  microwave  switches.  For 
example  an  output  voltage  of  0.010  from  Touchstone  to  represent  a  value  of  0.0103 
volts  can  cause  an  error  of  0.88  watts  of  diode  power  in  a  one  kilowatt  switch.  This 
number  may  look  small,  however  it  can  cause  an  error  of  30-40  C  in  diode  junction 
temperature  rise.  This  error  can  affect  the  decision  of  selection  of  a  proper  p-i-n  diode 
for  the  switch  design. 

D.       EVALUATION     OF     POWER    DISSIPATION     AND     HEATING    USING 
FORTRAN  PROGRAM 

To    avoid   some   disadvantages    of  the   use    of  the   Touchstone    nodal   voltage 

evaluation  for  power  computation,  we  may  search  for  another  CAD  program  or  create 

a  specific  computer  program  to  solve  this  problem.    Many  microwave  CAD  programs 

have  been  written  in  FORTRAN  IV  for  both  mainframe  and  minicomputers,  but  up  to 

now  there  is  no  program  written  to  solve  for  the  power  dissipation  in  a  microwave 

circuit  [Ref  7].    The  Fortran  language  requires  specific  formats  of  input  and  output 
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and  has  to  be  compiled  before  running  the  program.  Nonetheless,  the  advantages  to 
the  microwave  engineer  are  its  complex  number  notation,  and  double  precision 
capability.  These  two  advantages  lead  to  accurate  outputs,  for  example  in  the 
multiplication  of  more  than  ten  complex  matrices  in  cascade.  To  write  a  Fortran 
program  to  evaluate  the  terminal  voltages  of  the  diode  microwave  switch  equivalent 
circuit,  electrical  two-port  components  of  the  switch  are  converted  ^  into  simple 
appropriate  two-port  matrix  forms.  In  chapter  2  we  introduced  an  ABCD  matrix 
representation  of  the  two-port  networks  for  a  transmission  line  section,  a  shunt 
admittance  and  series  impedance  m  the  iumped-circuits.  For  our  switch  circuit  model, 
the  circuit  components  are  characterized  as  boxes  in  cascade  connecting  with  the 
applied  voltage  at  one  end  and  the  load  at  the  other  end. 

To  evaluate  the  power  dissipation  and  junction  temperature  in  the  diodes  of  a 
SPDT  switch  which  has  an  equivalent  circuit  shown  in  Figure  4.12,  with  the  optimal 
transmission  line  sections  by  use  of  Touchstone,  the  procedure  follows  the  flow 
diagram  shown  in  Figure  4.13. 

1.  Nodal  Voltage  and  Current  Calculation 

We  assume  that  the  switch  is  supplied  by  a  one  volt  generator  with  50  ohms 
series  impedance.  First  we  divide  the  equivalent  circuit  into  two  parts,  one  the  on- 
state  branch  the  other  the  off-state  branch.  For  both  on-  and  off-states,  each 
transmission  line  section  and  the  diode  equivalent  circuit  two-port  network  is 
transformed  into  ABCD  matrix  form  (Figure  4.12(a))  by  subroutine  MATRIX.  This 
subroutine  using  equations  2.18  and  2.19  for  the  transformation  of  a  shunt  diode  and  a 
lossless  transmission  line  section  into  ABCD  two-port  matrices  respectively.  To  obtain 
the  voltages  and  currents  in  the  circuit,  the  successive  input  impedances  of  the 
individual  two-port  elements  must  be  determined.  All  two-port  networks  in  the  branch 
are  cascaded  by  subroutine  MATMUL  and  produce  a  new  equivalent  circuit  as  in 
Figure  4.12(a).  Subroutine  MATMUL  provides  a  complex  matrix  product  of  a  pair  of 
matrices.  From  Figure  4.12(b),  input  impedance  in  each  branch  {'Z.-^^q^,  ^inofT^  ^^ 
calculated  from  equation  2.21.  The  impedance  of  the  switch  'Z-  )  is  determined  by  the 
paralleling  of  input  impeaances  of  both  branches  iZ-  _-, .  Z-^^r).  The  switch  input 
current  (I'^^)  and  input  voltage  (v-  )  from  the  generator  are  calculated  from  equations 
2.22  and  2.23.  The  next  step  is  to  find  the  currents  and  voltages  in  the  on-  and  off- 
state  branches.  The  input  voltage  on  each  branch  is  equal  to  the  switch  input  voltage 
(paralleled  ports).  Current  flow  in  the  on-  and  off-state  branches  is  determined  by 
current  divider  method  as  follows: 


60 


vvvwHi'' 


A^AA — fii 


N 


N 


r^ — ll^    r^==^'' 


> 


Figure  4.12    ABCD  Matrix  Representation  of  SPDT  switch  Circuit  Model. 
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SUBROUTINE  MATRIX 
ABCD  REPRZSZNTATION 


INPUT  PAKAiMETERS: 

1.  DIODE:  Rj,  Cj,  DIODE  Q,  NO. OF  OIODE 

2.  TRANSMISSION  LINE:  Zo,  EL 
:.  OPERATING  rRECUENC: 


Jt. 


SUBROUTINE  MATRIX 
ABCD  REPRESENTATION 


SUBROUTINE  MATMUL 
N-TOO  PORTS  CXSC\DE 


OFF^STATE   SRANOi 
INPUT   IMPEDANCE 


OFF-STATE  3RANCT 
INPUT  VOLTAGE 
INPUT  CURRENT 


INPUT  PO'.'fEa 
POWEn  IN  LOAD 


SUBROUTINE   POWER 

TERMINAL  VOLTAGE 

OIODE  POWER 

DISSIPATION 


CI   INPUT  SIGNAL 
SUBROUTI^fE  THERM  I 


NO 


NO 


SWITCi: 
INPUT  IMPEDANCE 
INPUT  VOLTAGE 
INPUT  CURRENT 


-i. 


SUBROUTINE  MATMUL 
N-TWO  PORTS  CASCADE 


ON-STATE  3RANCT 
INPUT  IMPEDANCE 


NO 


ON-STATE   BRANCi 
INPUT  VOLTAGE 
INPUT  CURRENT 


INPUT   POWER 
POlfES   IN  LOAD 


SUBROUTINE  POWER 

TERMINAL   VOLTAGE 

DIODE   POWER 

DISSIPATION 


PULSE  INPUT  SIGNALI 
UoRCUTI^fE  THE?-!- 


STOP 


Figure  4.13     Flow  Diagram  for  Calculation  of  Power  Dissipation 
and  Diode  Junction  Temperature. 
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I  =         ^in^^nofT  (eqn4.11) 


on 


^inon  "*"  ^inoff 


■^mon    ^inofr 

~  in "    on 

where   I^^  is  the  current  flow  in  switch  on-state  branch, 

\q^  is  the  current  flow  in  the  switch  ofi'-state  branch, 
^inon  ^^  ^^^  input  impedance  of  the  switch  on-state  branch  and, 
^inoff  ^^  ^^^  input  impedance  of  the  switch  off-state  branch. 
The  final  step  of  nodal  voltage  and  current  calculation  returns  to   Figure 
4.12(a).  In  each  branch,  starting  from  the  leftmost  node,  currents  and  voltages  in 
adjacent  nodes  are  calculated  by  equations  2.24  and  2.25.  In  this  step  we  obtain  the 
currents  flowing  in  and  out  of  each  particular  node,  as  well  as  the  input  and  output 
voltages,  leading  us  to  the  power  calculation  as  in  the  next  section. 

2.  Power  Dissipation  in  the  Individual  P-I-N  Diode 

For  the  calculation  of  diode  power  dissipation  in  each  branch  of  the  switch, 
we  create  the  subroutine  named  POWER.  In  this  subroutine,  the  input  currents  and 
input  voltage  of  each  branch  are  the  initial  input  parameters.  These  currents  together 
with  input  and  output  voltages  of  the  diode's  node  are  calculated  from  equations  2.24 
and  2.25.  Equation  2.26  is  used  to  calculate  the  diode  power  dissipation  with  its 
corresponding  nodal  voltage  and  current.  For  the  power  at  load  we  can  use  the 
rightmost  nodal  current  and  voltage  (I^^,  V^^)  as  inputs  to  equation  2.27  or  using 
equation  2. 28  with  the  branch  input  power  minus  the  total  of  all  power  dissipation  in 
that  branch  (include  power  dissipation  in  the  transmission  line  sections).  The  power 
dissipation  in  each  diode  and  all  the  calculated  powers  are  normalized  to  the  power 
delivered  to  the  load  in  the  switch  on-state. 

3.  Junction  Temperature  :n  the  Individual  P-I-N  Diode 

The  junction  temperature  of  a  p-i-n  diode  is  determined  by  the  ambient 
temperature  of  the  switch  and  the  power  dissipation  in  the  diode  itself  The  ambient 
temperature  is  assumed  to  be  constant  at  about  the  room  temperature.  The  power 
dissipation  in  the  diode  is  a  major  factor  that  causes  damage  to  the  semiconductor 
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devices  in  the  switch.  From  the  calculation  of  power  dissipation  mentioned  in  previous 

section,  we  can  separate  the  calculation  of  diode  junction  temperature  into  two  cases, 

CW  power  dissipation  and  pulsed  power  dissipation.  In  the  Fortran  program,  the 

subroutine  THERM!  evaluates  the  junction  temperature  of  the  CW-power  case,  for 

the  diode  mounted  on  a  heat  sink.  The  junction  temperature  is  considered  by  summing 

the  ambient  temperature  .and  the  increment  of  temperature  of  the  diode  caused  by  the 

input  CW-power  (equation  3.7).  The  total  thermal  resistance,  6:^,    of  a  particular 

diode  is  available  from  the  p-i-n  diode  specification  sheet.    For  the  pulsed  input-power 

case,  the  subroutine  THERM2  is  used  for  evaluating  the  diode  junction  temperature 

rise.  In  this  case  we  assume  that  the  input  pulse  width  is  small  compared  with  the 

thermal  time  constant,  meaning  that  all  the  thermal  energy  is  absorbed  in  the  diode, 

with   no   energy   spreading    out   to   the   surroundings.    The  junction   temperature   is 

obtained  by  summing  the  ambient  temperature  and  the  temperature  rise  of  the  junction 

(equation  3.12),  where  the  temperature  rise  of  the  junction  is  determine  by  the  thermal 

capacity  of  the  diode  active  region  (HC),  the  pulse  width  in  microsecond  and  the  power 

dissipation    in    the    diode    (equation    3.11).     For    the    switch    design    problem,    the 

recommended     maximum    junction     temperature     for     safe     operation     is     200  °C 

[Ref  10:  p.290]. 

4.  An  Example  and  Results  of  Power  Dissipation  and  Diode  Junction  Temperature 
Evaluation 

Appendix  D  is  the  listing  of  two  Fortran  programs  using  for  the  evaluation  of 

power  dissipation  and  junction  temperature  in  each  diode  of  a  single  pole  single  throw 

switch,  and  a  single  pole  double  throw  switch  with  one  or  more  diodes.    These  two 

programs,  PSPST  and  PSPDT,  are  written  for  general  use.  The  number  of  diodes  in 

each  branch  is  an  input  variable  that  determines  the  switch  construction  (number  of 

transmission  line  sections).  Diode  input  parameters  can  be  either  diode  Q  or  junction 

resistance  and  junction  capacitance.  All  input  parameters  can  be  entered  interactively 

or  in  an  unformatted  input  file,  and  the  results  of  the  evaluation  are  listed  in  an  output 

file.    As  an  example  of  power  dissipation  and  junction  temperature  in  a  SPDT. switch 

with  four  .\iA-4P203  p-i-n  diodes,  we  use  PSPDT  program  m  Appendix  D2  for  :he 

evaluation.    From   the  specificanon  sheet,   diode   MA-4P203   has  0.15   pF.  junction 

capacitance,  1.5  ohms  junction  resistance  and  30°C/watt  total  thermal  resistance.   The 

value    of    transmission    characteristic    impedance    and    electrical    length    for    each 

transmission    line     section    are    the     optimal    component    values     obtained    from 

SPDTD.CKT  Touchstone  circuit   file.     For  a    100  watt   CW  power   SPDT   switch 
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operating  at  10  to  20  GHz  frequency  range,  the  results  of  using  PSPDT  program  are 
listed  in  Table  9. 


TABLE  9 

DIODE  POWER  DISSIPATION  AND  JUNCTION  TEMPERATURE  OF 
THE  100- WATT  CW  POWER  SWITCH 

ON-STATE 

OFF-STATE 

PRE 

POWER(watts) 

TEMP(°C) 

POWER(watts) 

TEMP(°C) 

(GHz) 

Dl 

D2 

Dl 

D2 

Dl 

D2 

Dl 

D2 

10.00 

0.814 

0.578 

39.4 

32.3 

2.523 

0.002 

90.7 

15.1 

12.00 

1.124 

0.832 

48.7 

40.0 

2.653 

0.002 

94.6 

15.1 

14.00 

1.486 

1.133 

59.6 

48.9 

2.770 

0.002 

98.1 

15.1 

• 

16.00 

1.908 

1.479 

72.2 

59.4 

2.865 

0.002 

100.9 

15.1 

18.00 

2.407 

1.872 

87.2 

71.2 

2.932 

0.002 

102.9 

15.1 

20.00 

3.003 

2.311 

104.1 

84.3 

2.971 

0.002 

104.1 

15.1 

From  the  result  in  Table  9,  we  realize  that  most  of  the  switch  power  is 
dissipated  in  the  first  diode  of  each  branch,  particularly  in  the  switch-off  branch.  In  the 
switch  on- state  the  power  is  shared  among  the  two  diodes,  with  the  first  one  dissipating 
more  power.  For  the  whole  range  of  frequencies  the  first  diode  of  the  off-state  branch 
dominates  the  power  dissipation  and  undergoes  the  greatest  junction  temperature  rise. 
The  switch  can  be  damaged  if  the  junction  temperature  exceed  its  maximum  allowable 
value  (200  C  ).  In  most  diode  specification  sheets,  the  maximum  power  dissipation  at 
the  given  amoieni  temperature  is  given  based  on  the  CW-mput  power  to  the  diode. 
The  CW-input  power  has  more  effect  in  the  temperature  rise  in  the  diode  junction 
because  of  the  heat  accumulated  continuously  in  the  junction.  For  this  reason  a  diode 
switch  with  short-duration  input  power  pulses  can  handle  higher  peak  powers  than  for 
CW-input  power.  For  example  the  MA-4P203  diode  has  the  maximum  power  handling 
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capability  of  5  watts  of  CW-input  power  at  25  °C  ambient  temperature.  In  Table  9  the 
designed  switch  has  a  maximum  power  dissipation  in  the  diode  of  about  three  watts 
and  105''C  junction  temperature  (20*'C  ambient  temperature).  The  result  of  Table  9 
shows  that  our  switch  design  is  capable  to  handle  a  100-watt  input  CW  power.    Figure 

4.14  shows  the  plot  of  normalized  powers  to  one  watt  output  power  of  an  individual 
diode  in  both  branches. 
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POWER  DISSIPATION  VS  FREQUENCY 


12  14  13 

FREQUENCY  (GHZ) 


SPDT  4-DIODE  SWITCH 


Figure  4,14    Power  Dissipation  in  an  Individual  Diode. 
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V.  CONCLUSIONS  AND  REMARKS 

In  microwave  switching  design,  the  CAD  program  (Touchstone)  and  Fortran 
programs,  are  used  for  evaluating  different  types  of  microwave  diode  switches  in  order 
to  obtain  their  estimated  switching  characteristics.  The  selected  diode  Q  as  a 
switching-semiconductor  figure  of  merit  can  be  used  as  the  measure  of  device  switching 
effectiveness,  as  shown  in  chapter  4.  For  radar  applications,  a  single  pole  double  throw 
switch  with  high  isolation  over  a  broad  operating  frequency  range  is  an  essential 
requirement  for  all  radar  systems.  An  addition  of  more  diodes  to  the  switch  can 
improve  the  switch  isolation,  while  increasing  the  insertion  loss.  Figures  5.1  and  5.2 
are  plots  of  the  switch  insertion  loss  and  isolation  for  switches  with  one  or  more 
diodes.  The  insertion  loss  and  load  isolation  in  these  two  figures  are  produced  from  the 
Touchstone  circuit  file  evaluated  for  the  different  diode  Q  values.  These  two  figures 
are  very  useful  to  the  designer  selecting  the  number  of  diodes  to  be  used  for  specified 
switch  insertion  loss  and  isolation,  based  on  the  diode  Q  value.  An  appropriate  Q  value 
is  selected  to  characterize  the  behavior  of  this  switch.  The  Fortran  program  is  used  to 
estimate  the  power  handling  capability  of  the  designed  switch.  The  resistor  in  the  diode 
equivalent  circuit  plays  an  important  role  in  its  power  dissipation.  The  value  of  diode 
Q,  from  equation  4.1,  is  a  function  of  both  junction  capacitance  and  resistance.  One  of 
these  two  values  can  be  set  as  an  arbitrary  parameter,  meaning  that  we  can  select  more 
than  one  diode  having  the  same  diode  Q  value.  The  problem  is  to  select  that  p-i-n 
diode  which  can  handle  the  power  dissipation  of  the  switch.  One  suggestion  is  to  select 
lower  diode  resistance  to  avoid  too  much  power  dissipation.  However  this  leads  to  a 
trade-off  in  increased  power  dissipation  in  the  on-state  branch.  An  example  of  this  is  in 
the  previous  chapter;  the  p-i-n  diode  model  number  MA-4P203  has  the  junction 
resistance  of  1.5  ohms  and  diode  Q  of  47.15  at  15  GHz,  and  power  dissipation  in  the 
on-state  branch  of  the  fu-st  diode  is  very  close  to  that  of  the  first  diode  in  the  off-state 
branch  (1.68  and  2.82  watts).  In  some  switch  applications,  it  is  not  desirable  co  have 
heat  accumulated  in  the  diode  junction  at  all  times  (both  switching  state),  which  may 
cause  shorter  diode  lifetime.  If  we  select  a  p-i-n  diode  which  has  the  same  diode  Q 
value  but  greater  junction  resistance,  for  example  2.0  ohms  junction  resistance,  we 
obtain  3.478  watts  power  dissipated  in  the  off-state  branch  and  1.313  watts  in  the  on- 
state  branch.    The  best  way  to  select  the  diode  for  the  switch  is  by  using  the  Fortran 
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program  to  estimate  the  power  dissipation  in  the  first  diode  of  each  branch,  by  fixing 
the  diode  Q  value  and  varying  diode  junction  resistance  and  select  the  satisfactory 
outcome. 

Due  to  the  approximate  circuit  model  used  for  the  diode  microwave  switch 
equivalent  circuit,  the  power  calculation  cannot  produce  a  reliable  substitute  for 
measurements  in  the  laboratory.  Nonetheless  it  can  reduce  the  design  cost  by  providing 
initial  guide  lines  to  the  designer.  On  the  other  hand  these  assumptions  reduce  the 
complexities  of  the  circuit  analysis  and  are  mainly  employed  in  the  test  of  different 
figures  of  merit.    Improvements  for  more  accurate  result  can  be  achieved  as  follows: 

1.  Using  more  accurate  p-i-n  diode  equivalent  circuit  model. 

2.  Introducing  loss  in  the  transmission  line  sections. 

From  the  knowledge  that  the  first  diode  in  the  off-state  branch  dominates  the 
power  dissipation,  in  a  switch  circuit  design  for  more  efficient  switch  performance,  the 
diodes  used  in  the  circuit  can  be  differed  (not  necessarily  identical).  The  one  closest  to 
the  generator  should  have  higher  diode  Q  value  (lower  junction  resistance)  in  order  to 
dissipate  less  power,  allowing  the  remaining  diodes  to  share  the  power  dissipation. 
This  method  can  reduce  the  maximum  power  dissipation  that  occurs  in  the  first  diode 
of  the  off-state  branch,  and  improve  the  power  handling  capability  of  the  switch.  We 
still  can  use  the  Touchstone  circuit  file  to  evaluate  the  switch  performance  by  making 
some  necessary  modifications  to  it.  The  switch  optimal  insertion  loss  and  isolation  are 
processed  by  the  optimization  mode  to  obtain  the  optimal  switch  circuit  component 
values  to  be  used  in  the  power  dissipation  evaluation  with  Fortran  program. 
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Figure  5.1     Switch  Insertion  Loss  vs  Diode  Q  for  Multi-diodes  Switch. 
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Figure  5.2    Switch  Isolation  vs  Diode  Q  for  Mul'ti-diodes  Switch. 
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APPENDIX  A 
TOUCHSTONE  CIRCUIT  FILE  FOR  SPST  SWITCH 

This  appendix  provides  the  listing  of  two  circuit  PJes  using  in  the  Touchstone 
CAD  program  for  the  evaluation  switch  insertion  loss  and  load  isolation  for  single 
pole  single  throw  two-diode  switches. 

Al.  SPSTD.CKT  evaluate  with  diode  parameters  (C;,  R;). 
A2.  SPSTQ.CKT  evaluate  with  diode  figure  of  merit  (diode  Q). 

Al.  SPSTD.CKT 

SINGLE  POLE  SINGLE  THROW  DIODE  SWITCH 
INSERTION  LOSS  AND  LOAD  ISOLATION  CALCULATION 
INPUT  PARAMETERS  : 

1.  JUNCTION  CAPACITANCE,  CJ 

2.  JUNCTION  RESISTANCE,  RJ 

3.  CENTER  FREQUENCY,  FO 
OUTPUT  : 

OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE 

AND  ELECTRICAL  LENGTH 

MINIMIZED  INSERTION  LOSS 

LOAD  ISOLATION 
VAR 

!  INPUT  PARAMETERS: 

FO  =  15  !  FREQUENCY  IN  GHz 

CJ  =  .15  !  JUNCTION  CAPACITANCE  IN  pF 

RJ  =  1.5  !  JUNCTION  RESISTANCE  IN  Ohm 

!  OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE (Ohm) 
ZA#30  39.40  90 
ZB#30  38.25  90 

!  OPTIMIZED  TRANSMISSION  LINE  ELECTRICAL  LENGTH(Degree) 
EA#30  106.3  120 
EB#30  62.70  120 
CKT 
TLIN  1  2  Z'^ZA  E^EA  F^FO 
DEF2P  1  2  SEGl 
TLIN  1  2  Z'^ZB  E^EB  F'^FO 
DEF2P  1  2  SEG2 
SRC  1  2  R'^RJ  C^CJ 
DEF2P  1  2  DON 
RES  1  2  R'^RJ 
DEF2P  1  2  DOFF 
!  ON  STATE  CALCULATION 
SEGl  1  2 
DON  2  0 
SEG2  2  3 
DON  3  0 
SEGl  3  4 

DEF2P  1  4  LINDON 
!  OFF  STATE  CALCULATION 
SEGl  1  2 
DOFF  2  0 
SEG2  2  3 
DOFF  3  0 
SEGl  3  4 

DEF2P  1  4  LINDOFF 
OUT 
LINDON  DB(S21)  GRl  !  INSERTION  LOSS  IN  dB 
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LINDOFF  DB(S21)  GRIA 
FREQUENCY 

SWEEP  a  22  2 
GRID 

GRl   -2   0   .25 

GRIA  -20   -80 
OPT 

RANGE  10  20 

LINDCN  DB(321)=0 


!  LOAD  ISOLATION  IN  dB 

!  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 


!  MINIMIZE  INSERTION  LOSS 


A2.  SPSTQ.CKT 


SINGLE  POLE  SINGLE  THROW  DIODE  SWITCH 
INSERTION  LOSS  AND  LOAD  ISOLATION  CALCULATION 
INPUT  PARAMETERS  : 

1.  JUNCTION  CAPACITANCE,  CJ 

2.  FIGURE  OF  MERIT  (DIODE  Q) ,q 

3.  CENTER  FREQUENCY,  FO 
OUTPUT  : 

OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE 
AND  ELECTRICAL  LENGTH 
MINIMIZED  INSERTION  LOSS 
LOAD  ISOLATION 
VAR 

!  INPUT  PARAMETERS: 

FO  =  15  !  FREQUENCY  IN  GHz 

CJ  =  .15  !  JUNCTION  CAPACITANCE  IN  pF 

Q   =  100  !  FIGURE  OF  MERIT  (DIODE  Q) 

!  OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE (Ohm) 

ZA#30  36.66  90 

ZB#30  37.30  90 

!  OPTIMIZED  TRANSMISSION  LINE  ELECTRICAL  LENGTH (Degree) 

EA#30  108.3  120 

EB#30  58.7   120 
EQN 

RJ  =  1000/(2'^PI*F0*CJ*Q)      !  JUNCTION  RESISTANCE  CALCULATION 
CKT 

TLIN  1  2  Z'-^ZA  E'^EA  F^FO 

DEF2P  1  2  SEGl 

TLIN  1  2  Z'^ZB  E'^EB  F-^FO 

DEF2P  1  2  SEG2 

SRC  1  2  R'^RJ  C^CJ 

DEF2P  1  2  DON 

RES  1  2  R^RJ 

DEF2P  1  2  DOFF 

!  ON  STATE  CALCULATION 

SEGl  1  2 

DON  2  0 

SEG2  2  3 

DON  3  0 

SEGl  3  4 

DEF2P  1  4  LINDON 

!  OF?  STATE  CALCULATION 

SEGl  1  2 

DOFF  2  0 

SEG2  2  3 

DOFF  3  0 

SEGl  3  4 

DEF2P  1  4  LINDOFF 
OUT 

LINDON  DB(S21) 

LINDOFF  DB(S21) 
FREQUENCY 

SWEEP  8  22  2 
GRID 

GRl   -2  0   .25 


GRl 
GRIA 


!  INSERTION  LOSS  IN  dB 
!  LOAD  ISOLATION  IN  dB 

!  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 
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GRIA  -20   -80 
OPT 
RANGE  10  20 
LINDON  DB(S21)  =  0  !  MINIMIZE  INSERTION  LOSS 
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APPENDIX  B 
TOUCHSTONE  CIRCUIT  FILE  FOR  SPDT  SWITCH 

This  appendix  provides  the  listing  of  two  circuit  files  using  in  the  Touchstone 
CAD  program  for  the  evaluation  switch  insertion  loss  and  load  isolation  for  single 
pole  double  throw  four-diode  switches. 

Bl.  SPDTD.CKT  evaluate  \^ith  the  diode  parameters  (C:.  R:). 
B2.  SPDTQ.CKT  evaluate  with  diode  figure  of  merit  (diode  Q). 

Bl.  SPDTD.CKT 

! !  SINGLE  POLE  DOUBLE  THROW  DIODE  SWITCH 

! !  INSERTION  LOSS  AND  LOAD  ISOLATION  CALCULATION 

! !  INPUT  PARAMETERS  : 

!!  1.  JUNCTION  CAPACITANCE,  CJ 

!!  2.  JUNCTION  RESISTANCE,  RJ 

!!  3.  CENTER  FREQUENCY,  FO 

■ !  4 .  TERMINAL  LOAD 

!  OUTPUT  : 

!   OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE 

!   AND  ELECTRICAL  LENGTH 

!   MINIMIZED  INSERTION  LOSS 

!   LOAD  ISOLATION 
VAR 

!  INPUT  PARAMETERS 

FREQUENCY  IN  GHz 
JUNCTION  CAPACITANCE  IN  pF 
JUNCTION  RESISTANCE  IN  Ohm 
TERMINAL  LOAD  IN  Ohm 

!  OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE (Ohm) 

ZA#30  46.56  90 

ZB#30  39.06  90 

!  OPTIMIZED  TRANSMISSION  LINE  ELECTRICAL  LENGTH (Degree) 

EA#30  93.58  120 

EB#30  32.82  120 
CKT 

TLIN  1  2  Z^ZA  E^EA  F'^FO 

DEF2P  1  2  SEGl 

TLIN  1  2  Z'^ZB  E'^EB  F'^FO 

DEF2P  1  2  SEG2 

SRC  1  2  R'^RJ  C^CJ 

DEF2P  1  2  DON 

RES  1  2  R'^RJ 

DEF2P  1  2  DOFF 

!  ON  STATE  BRANCH 

SEGl  1  2 

DON  2  0 

SEG2  2  3 

DON  3  0 

SEGl  3  4 

!  OFF  STATE  BRANCH 

SEGl  1  5 

DOFF  5  0 

SEG2  5  6 

DOFF  6  0 

SEGl  6  7 

!  DEFINE  3  PORTS  SWITCH 

DEF3P  14  7  SWITCH 
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FO  =  15 
CJ  =  .15 
RJ  =  1.5 
RR  =   50 


SWITCH  12  3 

RES  3  0  R^RR 

DEF2P  1  2  SWON 

SWITCH  12  3 

RES  2  0  R-^RR 

DEF2P  1  3  SWOFF 
OUT 

SWON   DB(S21)  GRl 

SWOFF  DB(321)  GRIA 
FREQUENCY 

SWEEP  8  22  2 
GRID 

GRl   -2   0   .25 

GRIA  -20   -80 
OPT 

RANGE  10  20  ' 

LINDON  DB(S21)=0 


!  INSERTION  LOSS  SWITCH  ON  STATE 


!  LOAD  ISOLATION  SWITCH  OFF  STATE 

!  INSERTION  LOSS  IN  dB 

!  LOAD  ISOLATION  IN  dB 

!  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 


!  MINIMIZE  INSERTION  LOSS 


B2.  SPDTQ.CKT 


SINGLE  POLE  DOUBLE  THROW  DIODE  SWITCH 
INSERTION  LOSS  AND  LOAD  ISOLATION  CALCULATION 
INPUT  PARAMETERS  : 

1.  JUNCTION  CAPACITANCE,  CJ 

2.  FIGURE  OF  MERIT  (DIODE  Q) ,Q 

3.  CENTER  FREOUENCY,  FO 

4.  TERMINAL  LOAD,  RR 
OUTPUT  : 

OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE 
AND  ELECTRICAL  LENGTH 
MINIMIZED  INSERTION  LOSS 
LOAD  ISOLATION 


VAR 

INPUT  PARAMETERS 
FO  =  15 
CJ  =  .20 
Q   =  100 
RR  =  50 


FREQUENCY  IN  GHz 

JUNCTION  CAPACITANCE  IN  pF 

FIGURE  OF  MERIT  (DIODE  Q) 

TERMINAL  LOAD  IN  Ohm 
!  OPTIMIZED  TRANSMISSION  LINE  IMPEDANCE  (Ohm) 
ZA#30  36.66  90 
ZB#30  37.30  90 

!  OPTIMIZED  TRANSMISSION  LINE  ELECTRICAL  LENGTH  (Degree) 
EA#30  108.3  120 
EB#30  53.7   120 
EQN 

RJ  =  1000/(2*PI*F0'^CJ*Q)      !  JUNCTION  RESISTANCE  CALCULATION 
CKT 
TLIN  1  2  Z'^ZA  E^EA  F'^FO 
DEF2P  1  2  SEGl 
TLIN  1  2  Z^ZB  E'^EB  F'^FO 


DEF2P 
SRC  1 


2  SEG2 
R^RJ  C^CJ 


DHF2P  1  2  DON 

RES  1  2  R'^P.J 

DEF2?  1  2  DOFF 

!  ON  STATE  BRANCH 

SEGl  1  2 

DON  2  0 

SEG2  2  3 

DON  3  0 

SEGl  3  4 

!  OFF  STATE  BRANCH 

SEGl  1  5 

DOFF  5  0 

SEG2  5  6 
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DOFF  6  0 

SEGl  6  7 

!  DEFINE  3  PORTS  SWITCH 

DEF3P  14  7  SWITCH 

SWITCH  12  3 

RES  2  0  R-^RR 

DEF2P  1  2  SWON 

SWITCH  12  3 

RES  3  0  RARR 

DEF2P  1  3  SWOFF 
OUT 

SWON  DB(S21) 

SWOFF  DB(S21) 
FREQUENCY 

SWEEP  8  22  2 
GRID 

0   .25 
-80 


GRl 
GRIA 


GRl   -2 
GRIA  -20 
OPT 
RANGE  10 


20 


LINDON  DB(S21)  =  0 


!  INSERTION  LOSS  SWITCH  ON  STATE 


!  LOAD  ISOLATION  SWITCH  OFF  STATE 

!  INSERTION  LOSS  IN  dB 
I  LOAD  ISOLATION  IN  dB 

!  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 


!  MINIMIZE  INSERTION  LOSS 
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APPENDIX  C 
TOUCHSTONE  CIRCUIT  FILE  FOR  NODAL  VOLTAGE  EVALUATION 

This  appendix  provides  the  listing  of  two  circuit  files  using  in  the  Touchstone 
CAD  program  for  the  evaluation  of  nodal  voltages  by  using  voltage-controlled  voltage 
source  circuit  model  for  single  pole  single  throw  two-diode  switches  and  single  pole 
double  throw  four-diode  switches. 

CI.  SPSTV.CKT  single  pole  single  throw  nodal  voltage  evaluation 
C2.  SPDTV.CKT  single  pole  double  throw  nodal  voltage  evaluation 

CI.  SPSTD.CKT 

!  SINGLE  POLE  SINGLE  THROW  DIODE  SWITCH 

!  NODAL  VOLTAGE  CALCULATION  USING  VOLTAGE -CONTROLED  VOLTAGE  SOURCE 

!  INPUT  PARAI^IETERS  : 

!  1.  JUNCTION  CAPACITANCE,  CJ 

!  2.  JUNCTION  RESISTANCE,  RJ 

!  3.  CENTER  FREOUENCY ,  FO 

!  4.  TERMINAL  LOAD  IMPEDANCE,  RR 

!  5.  TRANSMISSION  LINE  IMPEDANCE,  ZA,  ZB 

!  6.  TRANSMISSION  LINE  ELECTRICAL  LENGTH,  EA,  EB 

!  OUTPUT: 

! !   NORMALIZED  NODAL  VOLTAGE  OF  THE  INTERESTED  NODE 
VAR 

!  INPUT  PARAMETERS: 

FO  =  15  !  FREQUENCY  IN  GHz 

CJ  =  .15  !  JUNCTION  CAPACITANCE  IN  pF 

RJ  =  1.5  !  JUNCTION  RESISTANCE   IN  Ohm 

RR  =  50  !  TERMINAL  LOAD  IMPEDANCE  IN  Ohm 

!  TRANSMISSION  LINE  IMPEDANCE (Ohm) 

ZA  =  36.66 

ZB  =  37.30 

!  TRANSMISSION  LINE  ELECTRICAL  LENGTH(Degree) 

EA  =  108.3 

EB  =  58.7 
CKT 

TLIN  1  2  Z'^ZA  E^EA  F^FO 

DEF2P  1  2  SEGl 

TLIN  1  2  Z^ZB  E^EB  F^FO 

DEF2P  1  2  SEG2 

CAP  1  2  C^CJ 

DEF2P  DCJ 

RES  1  2  R'^RJ 

DEF2P  1  2  DRJ 

!  ON  STATE  CALCULATION 

SEGl  1  2 

DCJ  2  5 

DRJ  5  0 

SEG2  2  3 

DCJ  3  6 

DRJ  6  0 

SEG  3  4 

RES  4  0  R-^RR 

VCVS  4  7  0  0  M=-.5  A=0  R1=1E8  R2=0  F=0  T=0 

VCVS  5  8  0  0  M=-.5  A=0  R1=1E8  R2=0  F=0  T=0 

VCVS  6  9  0  0  M=-.5  A=0  R1=1E8  R2=0  F=0  T=0 

DEF4P  17  8  9   OUTVl 
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!  OFF  STATE  CALCULATION 

SEGl  11  12 

DRJ  12  0 

SEG2  12  13 

DRJ  13  0 

SEGl  13  14 

RES  14  0  R'^RR 

VCVS  14  17  0  0  M=-.5 

VCVS  12  13  0  0  M=-.5 

VCVS  13  19  0  0  M=-.5 

DEF4P  11  17  18  19  0UTV2 
OUT 

!  NORMALIZED  NODAL  VOLTAGE 

OUTVl  MAG(S21)  GRl 

OUTVl  MAG(S31)  GRl a 

OUTVl  MAG(S41)  GRl A 

0UTV2  MAG(S2l)  GR2 

0UTV2  MAG(S3l)  GRIA 

0UTV2  MAG(S41)  GR2 
FREQUENCY 

SWEEP  8  22  2 
GRID 

RANGE  8   22   2 

GRl    0   .5   1 

GRIA  0   .15   .3 

GR2   0   .001   .0001 


A=0  R1=1E8  R2=0  F=0  T=0 

A=0  R1=1E8  R2=0  F=0  T=0 
A=0  Rl=iE8  R2=0  F=0  T=0 


!  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 


C2.  SPDTV.CKT 


! !  SINGLE  POLE  DOUBLE  THROW  DIODE  SWITCH 

! !  NODAL  VOLTAGE  CALCULATION  USING  VOLTAGE-CONTROLED  VOLTAGE  SOURCE 

! !  INPUT  PARAMETERS  : 

1.  JUNCTION  CAPACITANCE,  CJ 

2.  JUNCTION  RESISTANCE,  RJ 

3.  CENTER  FREQUENCY,  FO 

4.  TERMINAL  LOAD  IMPEDANCE,  RR 

5.  TRANSMISSION  LINE  IMPEDANCE,  ZA,  ZB 

6.  TRANSMISSION  LINE  ELECTRICAL  LENGTH,  EA,  EB 
!  OUTPUT  : 

!   NORMALIZED  NODAL  VOLTAGE  OF  THE  INTERESTED  NODE 
VAR 

!  INPUT  PARAMETERS: 

FO  =  15  !  FREQUENCY  IN  GHz 

CJ  =  .20  !  JUNCTION  CAPACITANCE  IN  pF 

RJ  =  1.5  !  JUNCTION  RESISTANCE  IN  Ohm 

RR  =  50  !  TERMINAL  LOAD  IMPEDANCE  IN  Ohm 

!  TRANSMISSION  LINE  IMPEDANCE (Ohm) 

ZA  =  35.90 

ZB   =  47.36 

!  TRANSMISSION  LINE  ELECTRICAL  LENGTH ( Degree ) 

HA   =  101.20 

Z3   =34.^ 

TLIN  1  2  Z^ZA  E'^EA  F'^FO 

DEF2P  1  2  SEGl 

TLIN  1  2  Z^ZB  E^EB  F'^FO 

DEF2P  1  2  SEG2 

CAP  1  2  C^CJ 

DEF2P  DCJ 

RES  1  2  R-^RJ 

DEF2P  1  2  DRJ 

!  ON  STATE  CALCULATION 

SEGl  1  2 
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DCJ  2  8 

DRJ  8  0 

SEG2  2  3 

DCJ  3  9 

DRJ  9  0 

SEG  3  4 

RES  4  0  R-^RR 

!  OFF  STATE  CALCULATION 

SEGl  1  5 

DRJ  5  0 

SEG2  5  6 

DRJ  6  0 

SEGl  6  7 

RES  7  0  R^RR 

VCVS  4  10  0  0  M=-. 5  A=0  R1=1E8  R2=0  F=0  T=0 

VCVS  8  11  0  0  M=-.5  A=0  R1=1E8  R2=0  F=0  T=0 

VCVS  5  12  0  0  M=-.5  A=0  R1=1E3  R2=0  F=0  1=0 

DEF4P  1  10  11  12  OUTVl 
OUT 

!  NORMALIZED  NODAL  VOLTAGE 

OUTVl  MAG(S21)  GRl 

OUTVl  MAG(S31)  GRl A 

OUTVl  MAG(S41)  GRl A 
FREQUENCY 

SWEEP  8  22  2  !  FREQUENCY  SWEEP  FROM  8  TO  22  GHz 

GRID 

RANGE  8   22   2 

GRl    0   .5   1 

GRIA  0   .5   .1 
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APPENDIX  D 

FORTRAN  PROGRAM  FOR  DIODE  POWER  AND  HEATING 

CALCULATIONS 

This  appendix  provides  tiie  listing  of  two  Fortran  programs  written  for  the 
evaluation  of  power  dissipation  and  junction  temperature  in  the  individual  diode  of  the 
single  pole  single  throw  and  single  pole  double  diode  switches. 

Dl.  PSPST  diode  power  dissipation  and  junction  temperature  of  a  single  pole 
single  throw  diode  switch. 

D2.  PSPDT  diode  power  dissipation  and  junction  temperature  of  a  single  pole 
double  throw  diode  switch. 

Dl.  PSPST 

CfiARACTER'*^!   ANSWER 

CHARACTER*13    FILE 

CHARACTER'*^3    SELECT 

CHARACTER'*=1    TYPE 

INTEGER  I,NPORT 

REAL  TA,TR,HC,TJON(10) ,TJOFF(10) ,DT 

REALMS  FRE , Q , C J , EL ( 10 ) , RR , PI , PIN , PON , POFF , PLON , PLOFF 

REALMS  PDON(IO) ,PDOFF(10) , AMP , VDONM,PEX 

C0MPLEX*16  VIN , ZOUT , Z ( 10 ) , ZINON , ZINOFF , ZIN , UN, ION , lOFF , ILON, VLON 

C0MPLEX*16  ILOFF,VLOFF,VON,VOFF,VS,ZS 

C0MPLEX*16  GPI(2,2,20),GP(2,2,20),GNI(2,2,20),GN(2,2,20) 

C0MPLEX*16  IDON ( 10 ) , IDOFF ( 10 ) , VDON ( 10 ) 

COMPLEX^ie  VDOFF(IO) 

OPEN (UNIT=2,FILE=' OUTPUT' ) 

PI  =  3.1415926536 
5      WRITE (6,*) 'POWER  DISSIPATION  AND  JUNCTION  TEMPERATURE  EVALUATION' 

WRITE(6,*) 'OF  SPST  DIODE  SWITCH  ' 

WRITE(6,*) '"PRESS  ANY  KEY  AND  RETURN"  TO  CONTINUE, "Q  AND  RETURN"TO 
*  QUIT' 

READ (5, 10) ANSWER 
10      FORMAT (Al) 

IF ( ANSWER. EQ. 'Q' )STOP 

WRITE(6,*) 'ENTER^'0"TO  READ  INPUT  FROM  FILE,  "1  "TO  ENTER  INPUTS  INTE 
*RACTIVELY' 

READ(5,*)NUM 

IF(NUM.NE.O.AND.NUM.NE.l)  GOTO  5 

IF(NUM.EQ.O)THEN 

MN  =  1 

ELSE 

MN  =  5 

SriDIF 

IF(NUM.£Q.O)THEN 

WRITE (6,*) 'ENTER  THE  NAME  OF  INPUT  FILE' 

READ(5,15)FILE 

0PEN(UNIT=1,FILE=FILE) 
15      FORMAT (Al 3) 

END  IF 

IF(NUM.EQ.1)THEN 
20     WRITE (6,*) 'WHAT  TYPE  OF  INPUT  SIGNAL  "PULSE"  OR  "CW"  ?  ' 

WRITE (6,*) 'FOR  "PULSE"  INPUT  ENTER  "P"  ' 

WRITE (6,*) 'FOR  "CW"  INPUT  ENTER  "C"  ' 

WRITE (6,*) 'DO  NOT  EVALUATE  JUNCTION  TEMPERATURE  ENTER  "N"  ' 
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.^?II?(§')!?S_™IS.  EVALUATION  USE  THE  FIGURE  OF  MERIT, DIODE  Q' 


ENDIF 

READ(MN,11)TYPE 
25     FORMAT(Al) 

IF(TYPE.NE. 'P' .AND. TYPE. NE. 'C .AND. TYPE .NE . 'N')GOTO  20 

IF(TYPE.EQ. 'N' )GOTO  30 

IF(TYPE.EQ. 'C )THEN 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  AMBIENT  TEMPERATURE  (DEGREE  C) ' 

READ(MN,*)TA 

IF(MUM.S0.1)WRITE(6,*) 'TOTAL  THERMAL  RESISTAI^ICS  (DEGREE  C/WATT)  ' 

READ(MN,^)TR 

ELSE 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  AMBIENT  TEMPERATURE  (DEGREE  C) ' 

READ(MN,'^)TA 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THERMAL  CAPACITY  OF  OF  I-REGION  (MICR 
*0-JOULE/DEGREE  C) ' 

READ(MN,''')HC 

IF (NUM. EQ.l) WRITE (6,^) 'ENTER  THE  PULSE  DURATION  (MICRO-SECOND)' 

READ(MN,^)DT 

ENDIF 
30      IF(NUM.EQ.1)THEN 

^;rite(6,*)  'IS  th: 

WRITE (6,*) 'ANSWER  "YES"  OR  "NO"' 
ENDIF 
READ (MN, 35) SELECT 
35     FORMAT (A3) 

IF (SELECT. NE. 'YES' . AND . SELECT .NE . 'NO' )GOTO  30 

IF(SELECT.EQ. 'YES' )THEN 

IF(NUM.EQ.iyWRITE(6,''f) 'ENTER  THE  VALUE  OF  DIODE  0  ' 

READ(MN,'')Q 

IF  (NUM.  EQ.l)  WRITE  (6,'^) 'ENTER  THE  VALUE  OF  JUNCTION  CAPACITANCE  (PF) 

READ(MN,'*=)CJ 

IF(NUM.E0.1)WRITE(6,*)'ENTER  THE  EVALUATE  FREQUENCY  (GHZ)' 

READ(MN,^)FRE 

RR  =  1000/(2'^PI*FRE*Q*CJ) 

ELSE 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  VALUE  OF  JUNCTION  CAPACITANCE  (PF) 

READ(MN,''^)CJ 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  EVALUATE  FREQUENCY  (GHZ)' 

READ(MN,^)FRE 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  VALUE  OF  JUNCTION  RESISTANCE  (OHM 

*)' 

READ(MN,*)RR 

Q  =  1000/(2*PI*FRE*RR'^CJ) 

ENDIF 

VS  =  (1,0) 

ZS  =  (50,0) 

IF (NUM. EO.l) WRITE (6,*) 'ENTER  NUMBER  OF  DIODES  IN  EACH  BRANCH' 

READ(MN,^)ND 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  EXPECT  POWER  AT  LOAD  (WATT)' 

READ(MN,'^)PEX 

IF(NUM.E0.1)WRITE(6,*) 'ENTER  THE  OUTPUT  IMPEDANCE  (COMPLEX)' 

READ(MN,'^)ZOUT 

NSEG  =  ND+1 

DO  40  I  =1,NSEG 

IFfNUM.E0.1)WRITS(6,*) 'TRANSMISSION  SEGMENT  NUMBER ', I 

IF(NUM.Zd.l;WRITE;6,^) 'ENTER  IMPEDANCE  OF  TRAI^ISMISSIOM  LINE  (COMPL 
^SX)  '    '         ' 

READ(MN,-^)Z(I) 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  ELECTRICAL  LENGTH  (DEGREE)' 

READ(MN,'^)EL(I) 

EL(I)  =  EL(I)*PI/180 
40      CONTINUE 

NPORT  =  ND+NSEG 

C      ON  STATE  CALCULATION 

CALL  MATRIX(RR,CJ,Z,EL, NSEG, ND,FRE,GP( 1,1,1)) 
DO  45  K  =1, NPORT 
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45 


IF(K.EQ.1)THEN 
GPI(1,1,1)  =  GP(1,1,1 
GPI( 1,2,1}  =  GP( 1,2,1 
GPI(2,l,l)  =  GP(2,1,1 
GPI(2,2,1)  =  GP(2,2,1 
ELSE 

CALL  MATMUL(GPI(1,1,K-1),GP(1,1,K),2,2,2,GPI(1,1,K)) 
END  IF 
CONTINUE 

ZINON  =  (GPI  (1 , 1  ,NPORT)'^ZOUT+GPI  (1 ,  2  ,NPORT)  )/  (GPI  (2  , 1  ,NPORT)'^ZOUT 
#+GPI(2,2,NPORT)) 


50 


OFF  STATE  CALCULATION 

CALL  MATRIX(RR,0.D0,Z,EL,NSEG,ND,FRE,GN(1,1,1)) 
DO  50  N  =1,  NPORT 
IF{N.EQ.1)THEN 
GNI( 1,1,1)  =  GN( 1,1,1' 
GNI(l,2,l)  =  GN(1,2,1' 
GNI(2,1,1)  =  GN(2,1,1' 
GNI(2,2,1)  =  GN(2,2,1 
ELSE 

CALL  MATMUL(GNI(1,1,N-1),GN(1,1,N),2,2,2,GNI(1,1,N)) 
ENDIF 
CONTINUE 

ZINOFF  = (GNI ( 1 , 1 , NPORT ) *ZOUT+GNI (1,2, NPORT ) ) / ( GNI ( 2 , 1 , NPORT ) ^ZOUT 
#+GNI(2,2,NP0RT)) 


ON-OFF   STATE 

VON      =      ZINON'^VS/(ZINON+VS) 

ION      =      VON/ZINON 

VOFF   =      ZINOFF*VS/(ZINOFF+VS) 

lOFF   =      VOFF/ZINOFF 

PON  =  0.5*REAL(VON'^DCONJG(ION)) 

POFF  =  0.5*REAL(VOFF*DCONJG(IOFF)) 

ILON  =  GPI  (1,1,  NPORT  )'^ION-GPI(  2,1,  NPORT  )*VON 

VLON  =  GPI(2,2,NPORT)'^VON-GPI(1,2,NPORT)*ION 

PLON  =  PON-0.5*REAL(VON'^DCONJG(  ION) -VLON*DCONJG(  ILON)) 

AMP   =  PEX/PLON 

WRITE (2, 55)  FRE 

WRITE(2,60)  Q 

WRITE (2, 65)  ND 

WRITE(2,70)  CJ 

WRITE (2, 75)  RR 

WRITE (2, 85)  PON*AMP 

WRITE (2, 90)  POFF^AMP 

WRITE (2, 95)  PLON* AMP 

ILOFF  =  GNI (1,1, NPORT )*IOFF-GNI( 2,1, NPORT )*VOFF 

VLOFF  =  GNI(2,2,NPORT)*VOFF-GNI(l,2,NPORT)'^IOFF 

PLOFF  =  POFF-0.5'^REAL(VOFF'^DCONJG(IOFF)-VLOFF*DCONJG(ILOFF)) 

WRITE (2, 100)  PLOFF^AMP 

'OPERATING  FREQUENCY' ,3X,F5. 2, 2X, 'GHZ' ) 

'FIGURE  OF  MERIT  (DIODE  Q)',3X,F8.3) 

'TOTAL  NUMBER  OF  DIODE  IN  EACH  BRANCH' , 2X, 12' 

' DIODE  JUNCTION  CAPACITANCE ' , 3X , F5 . 3 , 2X , ' PF ' 

' DIODE  JUNCTION  RESISTANCE ' , 3X , ?5 , 3 , 2X , ' OHM ' ) 

'POWER  IN  ON  STATE  ' , 2X , FIO . 2 , 3X, ' WATTS ' ^ 

'POWER  IN  OF?  STATE' , IX. ?10. 2, 3X, 'WATTS' ) 

'POWER  DELIVER  TO  THE  LOAD  IN  ON  STATE ', FIO . 2 , 3X, ' WATTS 

'POWER  DELIVER  TO  THE  LOAD  IN  OFF  STATE ', FIO .6 , 3X, ' WATT 


55 

FORMAT ( 

60 

FORMAT ( 

65 

FORMAT ( 

70 

FORMAT ( 

75 

FORMAT 

35 

FORMAT ( 

90 

FORMAT ( 

95 

FORI-IAT  { 
FORMAT ( 

100 

*S') 

POWER  DISSIPATED  IN  ON  STATE  BRANCH 

CALL  POWER ( GP ( 1 , 1 , 1 ) , ION , VON , NPORT , PDON , IDON , VDON ) 

WRITE (2,*) 'POWER  DISSIPATION  IN  DIODE  ON  STATE' 

DO  105  I  =  2,2*ND,2 

WRITE(2,115)I/2,PD0N(I)*AMP 
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105    CONTINUE 

C      POWER  DISSIPATED  IN  OFF  STATE  BRANCH 

CALL  POWER (GN( 1 ,1,1), lOFF , VOFF , NPORT , PDOFF , IDOFF , VDOFF) 

WRITE (2,*) 'POWER  DISSIPATION  IN  DIODE  OFF  STATE' 

DO  110  J  =  2,2*ND,2 

WRITE (2, 115) J/2, PDOFF (J)*AMP 
110     CONTINUE 

115    FORMATC  ', 'POWER  DISSIPATED  IN  DIODE  NO .  MX,  12  ,  2X,  FIO  .  5  ,  2X,  '  WATT 
'*'S'  ) 

C      JUNCTION  TEMPERATURE  CALCULATION 

C      JUNCTION  TEMPERATURE  IN  ON  STATE  BRANCH 

IF(TYPEoEQ. 'N' )STOP 

IF(TYPE.EQ. 'C )THEN 

CALL  THERMl (ND , T JON , PDON , AMP , TA , TR ) 

ELSE 

CALL  THERM2(ND,TJ0N,PD0N,AMP,TA,DT,HC) 

END  IF 

WRITE(2,*)' JUNCTION  TEMPERATURE  ON  STATE  BRANCH' 

DO  120  I  =  1,ND 

WRITE(2,82)I,TJ0N(I) 
120    CONTINUE 
C      JUNCTION  TEMPERATURE  IN  OFF  STATE  BRANCH 

WRITE (2,*) 'JUNCTION  TEMPERATURE  OFF  STATE  BRANCH' 

IF(TYPE.EO. 'C )THEN 

CALL  THERMl (ND,TJOFF, PDOFF, AMP, TA,TR) 

ELSE 

CALL  THERM2(ND,TJ0FF, PDOFF, AMP, TA,DT,HC) 

ENDIF 

DO  125  J  =  1,ND 

WRITE(2,130)J,TJOFF(J) 
125    CONTINUE 

130    FORMAT ('  ',' JUNCTION  TEMPERATURE  DIODE  NO. ' , 2X, 12 , 2X, F6 . 2 , 2X, ' DEGR 
*EE  C) 

STOP 

END 

C  SUBROUTINE  MATRIX 

SUBROUTINE  MATRIX (RR , C J , Z , EL , NSEG , ND , FRE , GP ) 

REALMS  RR,CJ,EL(10),FRE,PI 

C0MPLEX*16  Z{10),Z1 

COMPLEX^ie  G(2,2,20),G1(2,2,20),GP(2,2,20) 

PI  =  3.1415926536 

DO  10  K  =  1,NSEG 

G(1,1,K)  =  DCMPLX(COS(EL(K)),0.D0) 

G(1,2,K)  =  DCMPLX(0.D0,1.D0)*Z(K)*DCMPLX(SIN(EL(K)),0.D0 

G(2,1,K)  =  DCMPLX(0.D0,1.D0)/Z(K)*DCMPLX(SIN(EL(K)),0.D0 

G(2,2,K)  =  G(1,1,K) 
10      CONTINUE 

DO  20  L=1,ND 

IF(CJ.EQ.0.D0)THEN 

Zl  =  DCMPLX{RR,O.DO) 

ELSE 

Zl  =  DCMPLX(RR,-1000/(2*PI*FRE*CJ)) 

END  IF 

G1(1,1,L)  =  DC:-IPLX(1.D0,0.D0) 

G1(1,2,L)  =  DCMPLX( 0.00,0. DO) 

Gl(2,i,L)=  1/Zl 

G1(2,2,L)  =  G1(1,1,L) 
20      CONTINUE 

DO  30  M  =1,NSEG 

N  =  2*M-1 

GP(i,i,N)  =  g(i,i,m; 

GP(1,2,N)  =  G(1,2,M 
GP(2,1,N)  =  G(2,1,M 
GP(2,2,N)  =  G(2,2,m; 
30      CONTINUE 

DO  40  I  =1,ND 
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J  =  2*1 

GP(1,1,J)  =  G1(1,1,I' 
GP(1,2,J)  =  Gl  1,2,1 
GP(2,1,J)  =  G1(2,1,I 
GP(2,2,J)  =  G1(2,2,I' 
40     CONTINUE 
RETURN 
END 


C  SUBROUTINE  MATMUL 

SUBROUTINE  MATMUL ( AA, BB ,NN, MM, RR, CC) 

INTEGER  NN,MM,RR 

'COMPLEX^ie  AA(2,2,1),BB(2,2,1),CC(2,2,1) 

DO  10  I  =  1,NN 

DO  20  J  =  1,RR 

CC(I,J,1)  =  0 

DO  30  K  =1,MM 

CC(I,J,1)  =  CC(I,J,1)  +  AA(I,K,1)*BB(K,J,1) 
30      CONTINUE 
20     CONTINUE 
10      CONTINUE 

RETURN 

END 


C  SUBROUTINE  POWER 

SUBROUTINE  POWER(G, I , V,N,PD , IID, VID) 

REALMS  PD(IO) 

C0MPLEX*16  G(2,2,20),V,I,IID(20),VID(20),IOD(20),VOD(20) 

DO  10  J=2,N,2 

IID(J)  =  G(1,1,J-1)*I  -  G(2,1,J-1)*V 

VID(J)  =  G(2,2,J-1)*V  -  G(1,2,J-1)*I 

lOD(J)  =  G(1,1,J)*IID(J)  -  G(2,1,J)*VID(J) 

VOD(J)  =  G(2,2,J)*VID(J)  -  G(1,2,J)*IID(J) 

PD(J)   =  0.5*REAL(VID(J)*DCONJG(IID(J))-VOD(J)*DCONJG(IOD(J))) 

I  =  lOD'"^ 

V  =  VOD 
10      CONTINUE 

RETURN 

END 


I 


k 


\i] 


SUBROUTINE  THERMl 
SUBROUTINE  THERMl (ND,TJ ,PD , AMP ,TA,TR) 
INTEGER  ND 

REAL  TJ(10),AMP,TA,TR 
REALMS  PD(IO) 
DO  5  1=1, ND 

TJ(I)  =  TA  +  (PD(I*2)*AMP)*TR 
CONTINUE 
RETURN 
END 


SUBROUTINE  THERM2 
SUBROUTINE  THERM2 (ND , T J , ?D , AMP , TA , DT ,HC) 
INTEGER  ND 

REAL  TJ(10),AMP,TA,DT,HC 
REALMS  PD(IO) 
DO  5  J  =1,ND 

TJ(J)  =  TA  +  (PD(J*2)'^AMP)'*^DT/HC 
CONTINUE 
RETURN 
END 


85 


D2.PSPDT 


CHARACTER'^1    ANSWER 

CHARACTER*13    FILE 

CHARACTER'S    SELECT 

CHARACTER'*^!   TYPE 

INTEGER  I,NPORT 

REAL  TA,TR,HC,TJOM(10) ,TJOFF(10) ,DT 

REAL'S  FRE,0,CJ,EL(10) ,RR,PI ,PIN ,PON ,POFF ,PLON,PLOFF 

REAL'S  PDONUO)  ,PDOFF(10)  ,  AMP  ,  VDONM,  PEX 

COMPLEX'16  VIN , ZOUT , Z ( 10 ) , ZINON , ZINOFF , ZIN , UN , ION , lOFF , ILON , VLON 

COMPLSX'16  ILOFF,VLOFF,VON,VOFF,VS,ZS 

COMPLEX'16  GPI ( 2 , 2 , 20 ) , GP ( 2 , 2 , 20 ) , GNI ( 2 , 2 , 20 ) , GN ( 2 , 2 , 20 ) 

COMPLEX'16  IDON ( 10 ) , IDOFF ( 10 ) , VDON ( 10 ) 

COMPLEX'16  VDOFF(IO) 

0PEN(UNIT=2,FILE=' OUTPUT' ) 

PI  =  3.1415926536 
5      WRITE (6,') 'POWER  DISSIPATION  AND  JUNCTION  TEMPERATURE  EVALUATION' 

WRITE(6,') 'OF  SPDT  DIODE  SWITCH  ' 

WRITE(6,') '"PRESS  ANY  KEY  AND  RETURN"  TO  CONTINUE, "Q  AND  RETURN"TO 
'  QUIT' 

READ (5, 10) ANSWER 
10     FORMAT (Al) 

IF(ANSWER.EQ. 'Q' )STOP 

WRITE(6,') 'cNTER"0"TO  READ  INPUT  FROM  FILE,"1"T0  ENTER  INPUTS  INTE 
'RACTIVELY' 

READ(5,')NUM 

IF(NUM.NE.O.AND.NUM.NE.l)  GOTO  5 

IF(NUM.EQ.O)THEN 

MN  =  1 

ELSE 

MN  =  5 

END  IF 

IF(NUM.EQ.O)THEN 

WRITE (6,') 'ENTER  THE  NAME  OF  INPUT  FILE' 

READ(5,15)FILE 

0PEN(UNIT=1,FILE=FILE) 
15     F0RMAT(A13) 

END  IF 

IF(NUM.EQ.1)THEN 
20     WRITE (6 ,')' WHAT  TYPE  OF  INPUT  SIGNAL  "PULSE"  OR  "CW"  ?  ' 

WRITE (6,') 'FOR  "PULSE"  INPUT  ENTER  "P"  ' 

WRITE (6,') 'FOR  "CW"  INPUT  ENTER  "C"  ' 

WRITE (6,') 'DO  NOT  EVALUATE  JUNCTION  TEMPERATURE  ENTER  "N"  ' 

END  IF 

RE AD (MN, 11) TYPE 
25      FORMAT (Al) 

IF(TYPE.NE. 'P' .AND. TYPE. NE. 'C . AND . TYPE . NE . 'N' )GOTO  20 

IF (TYPE. EO. 'N' )GOTO  30 

iF(TYPE.EO.  'c  When 

I? (NUM. HO Ti) WRITE (6,^) 'ENTER  AMBIENT  TEMPERATURE  ^DEGREE  C) ' 
READ(!-m,^)TA 

IF(NUM.E0.i)WRITE(6,^) 'TOTAL  THERMAL  RESISTANCE  (DEGREE  C/WATT) ' 
READ(MN,')TR 
ELSE 

IF(NUM.EQ.1)WRITE(6,')'ENTER  AMBIENT  TEMPERATURE  (DEGREE  C) ' 
READ(MN,')TA 

IF(NUM.EQ.1)WRITE(6,') 'ENTER  THERMAL  CAPACITY  OF  OF  I-REGION  (MICR 
'0-JOULE/DEGREE  C) ' 
READ(MN,')HC 

IF(NUM.EQ.1)WRITE(6,') 'ENTER  THE  PULSE  DURATION  (MICRO-SECOND)' 
READ(MN,'^)DT 
END  IF 
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30      IF(NUM.EQ.1)THEN 

WRITE(6,*) 'IS  THIS  EVALUATION  USE  THE  FIGURE  OF  MERIT, DIODE  Q 
^.  ^y -_  


WRITE (6,*) 'ANSWER  "YES"  OR  "NO" 
END  IF 
READ (MN, 35) SELECT 
3  5     FORMAT ( A3 ) 

IF ( SELECT. NE. 'YES' .AND. SELECT. NE. 'NO') GOTO  30 

IF(SELECT.EO. 'YES' )THEN 

IF(NUM.S0.lTWRITE(6,'^) 'ENTER  THE  VALUE  OF  DIODE  0  ' 

READ(MN,^)Q 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  VALUE  OF  JUNCTION  CAPACITAI^ICE  (PF) 

READ(MN,'^)CJ 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  EVALUATE  FREQUENCY  (GHZ)' 

READ(MN,*)FRE 

RR  =  1000/(2^PI*FRE*0*CJ) 

ELSE 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  VALUE  OF  JUNCTION  CAPACITANCE  (PF) 

READ(MN,*)CJ 

IF(NUM'.EQ.1)WRITE(6,*) 'ENTER  THE  EVALUATE  FREQUENCY  (GHZ)' 

READ(MN,'^)FRE 

IF(NUM.EQ.1)WRITE(6,=^) 'ENTER  THE  VALUE  OF  JUNCTION  RESISTANCE  (OHM 

*)' 

READ(MN,''^)RR 

Q  =  1000/(2*PI*FRE'^RR*CJ) 

END  IF 

VS  =  (1,0)^ 

ZS  =  (50,0) 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  NUMBER  OF  DIODES  IN  EACH  BRANCH' 

READ(MN,^)MD 

IF(NUM.EQ.1)WRITE(6,'^) 'ENTER  THE  EXPECT  POWER  AT  LOAD  (WATT)' 

READ(MN,^)PEX 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  THE  OUTPUT  IMPEDANCE  (COMPLEX)' 

READ(MN,'^)ZOUT 

NSEG  =  ND+1 

DO  40  I  =1,NSEG 

IF(NUM.EQ.1)WRITE(6,*) 'TRANSMISSION  SEGMENT  NUMBER', I 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  IMPEDANCE  OF  TRANSMISSION  LINE  (COMPL 
'^EX)  ' 

READ(MN,*)Z(I) 

IF(NUM.EQ.1)WRITE(6,*) 'ENTER  ELECTRICAL  LENGTH  (DEGREE)' 

READ(MN,^)EL(I) 

EL(I)  =  EL(I)*PI/180 
40      CONTINUE 

NPORT  =  ND+NSEG 

C      ON  STATE  CALCULATION 

CALL  MATRIX ( RR , C J , Z , EL , NSEG , ND , FRE , GP ( 1 , 1 , 1 ) ) 

DO  45  K  =1, NPORT 

IF(K.EQ.1)THEN 

GPI(1,1,1)  =  GP(1,1,1 

GPI(1,2,1)  =  GP(1,2,1 

GPI(2,1,1)  =  GP(2,1,1 

GPI(2,2,1)  =  GP(2,2,1' 

ELSE 

CALL  MATMUL(GPI(1,1,K-1) , G?(l , 1 ,K) , 2 , 2 , 2 , GPI( 1 , 1 ,K) ) 

ENDIF 
45      CONTINUE 

ZINON  =  (GPI  (1,1, NPORT) '^ZOUT+GPI(  1,2, NPORT)  )/(GPI(  2,1, NPORT) ''^ZOUT 
#+GPI(2,2,NP0RT)) 

C      OFF  STATE  CALCULATION 

CALL  MATRIX ( RR , 0 . DO , Z , EL , NSEG , ND , FRE , GN ( 1 , 1 , 1 ) ) 

DO  50  N  =1,  NPORT 

IF(N.EQ.1)THEN 

GNI(1,1,1   =  GN(1,1,1) 

GNI(1,2,1)  =  GN(1,2,1) 
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50 


GNI(2,1,1)  =  GN(2,1,1) 
GNI(2,2,1)  =  GN(2,2,1) 
ELSE 

CALL  MATMUL(GNI(1,1,N-1),GN(1,1,N),2,2,2,GNI(1,1,N)) 
END  IF 
CONTINUE 

ZINOFF  =(GNI(l,l,NPORT)*ZOUT+GNI(l,2,NPORT))/(GNI(2,l,NPORT)*ZOUT 
#+GNI(2,2,NPGRT)) 
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60 
65 
70 
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80 
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ON-OFF 

ZIN  = 

VIN  = 

UN  = 

PIN  = 

ION  = 

lOFF  = 

PON   = 

POFF  = 

ILON  = 

VLON  = 

PLON  = 

AMP   = 

WRITE  " 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

WRITE 

ILOFF 

VLOFF 

PLOFF 

WRITE ( 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 

FORMAT 
*S') 


STATE 
ZINON*ZINOFF/ (ZINON+ZINOFF) 
ZIN'^VS/(ZIN+ZS) 
VIN/ZIN 

REAL ( . 5* ( VIN^DCONJG ( I IN) ) ) 
IIM^ZINOFF/ ^ZINON+ZINOFF) 
UN  -  ION 

0 . 5*REAL ( VIN^DCON JG ( ION ) ) 
0.  S'^REAL  ( VIN'^DCON  JG  ( lOFF ) ) 
GPI  ( 1 , 1 , NPORT)  *ION-GPI  (2,1, NPORT) '^^VIN 
GPI (2,2, NPORT ) *VIN-GPI (1,2, NPORT ) *ION 
PON-0.5*REAL(VIN*DCONJG(ION)-VLON*DCONJG(ILON)) 
PEX/PLON 
2,55)  FRE 

Q 

ND 
CJ 
RR 

PIN^AMP 
PON^AMP 
^OF^'^-.MP 
PLON'^AMP 
GNI(1,1,NP0RT)*I0FF-GNI(2,1,NP0RT)*VIN 
GNI ( 2 , 2 , NPORT ) *VIN  -GNI ( 1 , 2 , NPORT ) *IOFF 
POFF-0  .  5*REAL  (VIN'^DCONJG  ( lOFF )  -VLOFF*DCONJG(  ILOFF)  ) 
100)  PLOFF^AMP 

'OPERATING  FREQUENCY ' , 3X , F5 . 2 , 2X , ' GHZ ' ) 
'FIGURE  OF  MERIT  (DIODE  Q)',3X,F8.3) 
'TOTAL  NUMBER  OF  DIODE  IN  EACH  BRANCH ', 2X, 12) 
JUNCTION  CAPACITANCE ' , 3X , F5 . 3 , 2X , ' PF ' ) 
JUNCTION  RESISTANCE' ,3X.F5.3,2X, 'OHM' ) 


2,60^ 

2,65 

2,70 

2,75 

2,80 

2,85 

2,90 

2,95) 


'DIODE 
'DIODE 

'INPUT  POWER 
'POWER  IN  ON 


'POWER 
' POWER 


,F10.2,3X, 'WATTS'  ] 

STATE  BRANCH ',2X,F10.2,3X,' WATTS ' ) 
IN  OFF  STATE  BRANCH ', IX, FIO . 2 , 3X, ' WATTS ' ) 
DELIVER  TO  THE  LOAD  IN  ON  STATE ', FIO . 2 , 3X, ' WATTS 


('  ',' POWER  DELIVER  TO  THE  LOAD  IN  OFF  STATE ', FIO . 6 , 3X, ' WATT 


C      POWER  DISSIPATED  IN  ON  STATE  BRANCH 

VON  =  VIN 

CALL  POWER ( GP ( 1 , 1 , 1 ) , ION , VON , NPORT , PDON , IDON , VDON) 

WRITE  (2,'^) 'POWER  DISSIPATION  IN  DIODE  ON  STATE' 

DO  105  I  =  2,2*ND,2 

WRITE(2,115)I/2,?D0N(I)*AI'1P 
105    CONTINUE 

C      POWER  DISSIPATED  IN  OFF  STATE  BRAI^ICH 

VOFF  =  VIN 

CALL  POWER (GN (1,1,1), lOFF , VOFF , NPORT , PDOFF , IDOFF , VDOFF) 

WRITE  (2,'*^) 'POWER  DISSIPATION  IN  DIODE  OFF  STATE' 

DO  110  J  =  2,2^ND,2 

WRITE (2, 115) J/2, PDOFF ( J) *AMP 
110     CONTINUE 

115     FORMAT('  ', 'POWER  DISSIPATED  IN  DIODE  NO.  MX,  12, 2X, FIO  .  5  ,  2X, 'WATT 
*S'  ) 

C      JUNCTION  TEMPERATURE  CALCULATION 

C      JUNCTION  TEMPERATURE  IN  ON  STATE  BRANCH 
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IF (TYPE. EQ. 'N' )STOP 

IF (TYPE. EQ. 'C )THEN 

CALL  THERMl (ND , T JON , PDON , AMP , TA , TR) 

ELSE 

CALL  THERM2(ND,T JON, PDON, AMP, TA,DT,HC) 

END  IF 

WRITE (2,*) 'JUNCTION  TEMPERATURE  ON  STATE  BRANCH' 

DO  120  I  =  1,ND 

WRITEf2,32)I,TJON(l) 
120     CONTINUE 
C      JUNCTION  TEMPERATURE  IN  OFF  STATE  BRANCH 

WRITE (2,*) 'JUNCTION  TEMPERATURE  OFF  STATE  BRANCH' 

IF (TYPE. EQ. 'C')THEN 

CALL  THERMl (ND , T JOFF , PDOFF , AMP , TA , TR) 

ELSE 

CALL  THERM2 (ND , T JOFF , PDOFF , AMP , TA , DT , HC ) 

END  IF 

DO  125  J  =  1,ND 

WRITE(2,130)J,TJOFF(J) 
125    CONTINUE 

130    FORMAT ('  ',' JUNCTION  TEMPERATURE  DIODE  NO . ' , 2X, 12 , 2X,F6 . 2 , 2X, ' DEGR 
'^EE  C  ) 

STOP 

END 


c  SUBROUTINE  MATRIX 

SUBROUTINE  MATRIX ( RR , C J , Z , EL , NSEG , ND , FRE , GP ) 

REAL  =^8  RR,CJ,EL(10)  ,FRE,PI 

COMPLEX^ie  Z(10),Z1 

C0MPLEX*16  G(2,2,20),G1(2,2,20),GP(2,2,20) 

PI  =  3.1415926536 

DO  10  K  =  1,NSEG 

G(1,1,K)  =  DCMPLX(COS(EL(K)),0.D0) 

G(1,2,K)  =  DCMPLX(0.D0,1.DG)'''Z(K)^DCMPLX(SIN(EL(K)),0.D0) 

G(2,1,K)  =  DCMPLX(0.D0,1.D0)/Z(K)'^DCMPLX(SIN(EL(K)),0.D0) 

G(2,2,K)  =  G(1,1,K) 
10      CONTINUE 

DO  20  L=1,ND 

IF(CJ.EQ.O.DO)THEN 

Zl  =  DCMPLX(RR,O.DO) 

ELSE 

Zl  =  DCMPLX(RR,-1000/(2*PI'^FRE*CJ)) 

ENDIF 

G1(1,1,L)  =  DCMPLX(1.DG,0.D0) 

G1(1,2,L)  =  DCM?LX(O.DO,O.DO) 

G1(2,1,L)=  1/Zl 

G1(2,2,L)  =  G1(1,1,L) 
20      CONTINUE 

DO  30  M  =1,NSEG 

N  =  2*M-1 

GP(i,i,N)  =  g(i,i,m; 

GP(1,2,N)  =  G(1,2,M 

GP(2,1,N)  =  G(2,1,M 

GP(2,2,N)  =  G(2,2,M 
30      CONTINUE 

DO  40  I  =1,ND 

J  =  2^1 

GP(1,1,J)  =  G1(1,1,I 

GP(1,2,J)  =  G1(1,2,I 

GP(2,1,J'^  =  G1(2,1,I 

GP(2,2,J)  =  G1(2,2,I 
40      CONTINUE 

RETURN 

END 
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C  SUBROUTINE  MATMUL 

SUBROUTINE  MATMUL (AA,BB,NN, MM, RR, CC) 

INTEGER  NN,MM,RR 

C0MPLEX*16  AA(2,2,1),BB(2,2,1),CC(2,2,1) 

DO  10  I  =  1,NN 

DO  20  J  =  1,RR 

CC(I,J,1)  =  0 

DO  30  K  =1,MM 

CC(I.J,1)  =  CCa,J,l)  +  AA(I,K,l)'^BB(K,J,l) 
30     CONTINUE 
20     CONTINUE 
10      CONTINUE  . 

RETURN 

END 


C  SUBROUTINE  POWER 

SUBROUTINE  POWER(G, I , V,N,PD , IID , VID) 

REALMS  PD(IO) 

C0MPLEX^16  G(2,2,20),V,I,IID(20),VID(20),IOD(20),VOD(20) 

DO  10  J=2,N,2 

IID(J)  =  G(1,1,J-1)*I  -  G(2,1,J-1)*V 

VID(J)  =  G(2,2,J-1)'^V  -  G(l,2,J-lpI 

lOD(J)  =  G(1,1,J)*IID(J)  -  G(2,1,J)*VID(J) 

VOD(J)  =  G(2,2,J)*VID(J)  -  G(1,2,J)'^IID{J) 

PD(J)   =  0.5^REAL(VID(J)*DCONJG(IID(J))-VOD(J)*DCONJG(IOD(J))) 

I  =  lOD(J) 

V  =  VOD(J) 
10      CONTINUE 

RETURN 

END 


SUBROUTINE  THERMl 
SUBROUTINE  THERI^l  (ND  ,  T J  ,  PD  ,  AMP ,  TA ,  TR ) 
INTEGER  ND 

REAL  TJ(IO) ,AMP,TA,TR 
REAL*8  PD(IO) 
DO  5  1=1, ND 

TJ(I)  =  TA  +  (PD(I*2)*AMP)'^TR 
CONTINUE 
RETURN 
END 


SUBROUTINE  THERM2 
SUBROUTINE  THERM2 (ND , T J , PD , AMP , TA , DT , HC ) 
INTEGER  ND 

REAL  TJ(IO) ,AMP,TA,DT,HC 
REALMS  PD(IO) 
DO  5  J  =1,ND 

TJ(J)  =  TA  +  (PD(J*2)*AMP)*DT/HC 
CONTINUE 
RETURN 
END 
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